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ABSTRACT 
The  transient  testing  technique  for  the  determination  of  heat  trans- 
fer characteristics  of  a  porous  medium  has  been  known  since  1929.   It 
requires  a  good  deal  of  care  and  one  of  the  important  restrictions  is  a 
step  change  in  temperature  must  be  applied  to  the  medium  at  the  beginning 
of  a  test  run.   This  thesis  involves  the  design,  test  and  evaluation  of 
an  experimental  set-up  that  utilizes  an  "integral  heater"  to  accomplish 
the  step  change  in  temperature.   In  addition,  there  is  presented  heat 
transfer  and  flow  friction  characteristics  for  a  20°  skew  matrix  con- 
structed from  perforated  nickel  material.   This  data  is  compared  to  a 
20°  skew  matrix  made  from  stainless  steel  material  that  has  been  pre- 
viously reported  by  C.  P.  Howard.   Also,  data  is  presented  for  two  "paral- 
lel plate"  matrices,  one  constructed  from  brass  sheet,  the  other  from 
perforated  nickel  material.   A  comparison  of  their  relative  performances 
is  made. 
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NOMENCLATURE 

2 
Matrix  heat  transfer  area,  ft 

2 
Flew  cross-section  area,  ft 

2 
Matrix  frontal  area,  ft 


Solid  matrix  cross- section  area  available  for  thermal 
conduction,  ft 

2 
Conduction  area  corrected  for  effect  of  slots  ft 

Specific  heat  of  fluid  at  constant  pressure,  BTU/ib°F 

Matrix  material  specific  heat,  BTu7lb°F 

Fluid  heat  capacitance,  mc  ,  BTU/*F  sec 

P 

W 

s 

E  ,   =    Hydraulic  Diameter,  4r  ,  ft 

£    =    Flow  friction  power  per  unit  area,  hp/ft 

2 
G    =    Mass  flow  velocity,  ib/hr  ft 

2 
g    -    Proportionality  constant,  32.20  ft  lb/#  sec 

h    =    Unit  conductance  for  connective  heat  transfer  or  heat  transfer 
power  per  unit  area  per  degree  temperature  difference,  BTU/hr 
CF  ft2 

k  ■  Fluid  thermal  conductivity,  BTU/hr°F  ft 

k  =  Matrix  material  thermal  conductivity,  BTU/hr°F  ft 

L  =  Matrix  Flow  length,  ft 

a  =  Mass  flow  rate,  Ib/hr 

V  :  pressure,  Ib/hr 

r    =    Hydraulic  radius,  ft 

V  •     Matrix  Volume  ft 
m 

V  -    Material  volume  corrected  for  effect  of  slots 
s 

3 

V  -         Material  volume  ft 


W    =    Matrix  mass,  lb 
s 

jS  -  Matrix  compactne 

jS  -  Compactness  corrected  for  effect  of  slots 

JUL  -  Fluid  viscosity,  lb/hr  ft 

a  s  Heat  transfer  area  reduction  factor 

q       -         Conduction  area  reduction  factor 

3 
P      ■    Fluid  density,  lb/ft 

0  m         Matrix  density,  lb/ft3 

Nondimensional  parameters: 

f    a    Fanning  friction  factor,  ratio  of  the  wall  shear  stress 
to  the  fluid  dynamic  head 

j    =    Colburn  j  factor,  NgTNpR2/3 

N„„  b    Nusselt's  number,  h  D„, 
Nu  H 

k 
N    b    Reynold's  number  4r  G/ 

N.„,  b    Stanton  number,  h/Gc 
ST  p 

N    =    Prandtl  number,  /Ac   /k 

N^   =    Number  of  heat  transfer  units,  hA/mc 
tu  p 

K    =    Entrance  flow  coefficient  for  the  matrix 
e 

K    b    Exit  flow  coefficient  for  the  matrix 
c 

p    =    porosity  of  the  matrix 

Matrix  flow  void  volume/Matrix  volume 

7^   =    Conduction  parameter,  k  A  /Lmc 
'   s  s    p 

Subscripts : 

i    =    initial 

1  =    upstream 

2  =    downstream 
f    =    fluid 

viii 


s    a    solid  (matrix  material) 

STD  =    Standard  temperature  and  pressure 


1,    INTRODUCTION. 

A  compact  heat  exchanger  is  one  not  necessarily  of  small  bulk  and 
weight,  but  one  which  incorporates  a  heat  transfer  surface  having  a  high 
area  density. 

Arbitrarily,  it  has  been  specified  that  JO    —   200  ft  /ft  makes 

for  a  compact  heat  exchanger,  where  p     is  the  ratio  of  the  heat  trans- 

2    3      + 
fer  area  to  core  volume,  ft  /ft  .  (15) 

Motivation  for  using  compact  surfaces  is  to  gain  high  exchanger 
effectiveness  in  given  box  volume  and  weight  limitations  and  at  a  reason- 
able cost.   Great  impetus  has  been  given  the  compact  heat  exchanger  field 
by  recent  applications  as  a  regenerator  in  aircraft  turboprop  power  plants, 
as  a  regenerator  in  the  gas  turbine  automobile  engine  and  by  their  use 
in  cryogenic  plant  exchangers  in  the  production  of  liquid  oxygen  and 
nitrogen.   As  the  applications  for  compact  heat  exchangers  increase, 
design  problems  increase.   A  successful  compact  heat  exchanger  design 
and  development  program  requires  access  to,  or  the  ability  to  generate 
knowledge  of  the  basic  heat  transfer  and  friction  behavior  of  compact 
surfaces  of  interest.   It  is  also  desirable  that  the  designer  would  have 
available  as  much  information  on  surface  performances  as  possible  to  re- 
main flexible  in  his  design  with  regard  to  the  limitations  that  have  been 
placed  upon  him. 

Building  a  complete  heat  exchanger  to  evaluate  a  new  surface  as  to 
its  heat  transfer  and  friction  performance  would  be  prohibitively  ex- 
pensive.  The  alternate  method  is  to  generate  the  performance  character- 
istics of  a  surface  of  interest  by  testing  small  matrix  modules  and  use 

+Numbers  in  parentheses  refer  to  references  in  the  bibliography. 
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this  data  to  design  a  prototype  heat  exchanger. 

Briefly,  the  objectives  of  this  thesis  are:   to  build  an  experimental 
rig  that  will  test  small  matrix  modules  to  determine  their  heat  transfer 
and  flow  friction  characteristics,  and,  to  test  matrices  made  from  per- 
forated nickel  material  to  determine  these  characteristics. 

One  of  the  simplest  and  most  satisfactory  methods  of  measuring  the 
forced  convection  heat  transfer  characteristics  of  a  matrix  utilizes  a 
transient  test  technique  or  so  called  single  blow  method.   This  test 
technique  consists  of  heating  the  homogeneous  porous  solid  of  known 
weight  and  geometry  with  a  uniformly  distributed  heat  source  and  then 
at  a  given  instant,  removing  the  heat  source  and  simultaneously  cooling 
the  solid  with  a  fluid  at  a  uniform  lower  temperature.   The  heating  and 
cooling  steps  can  be  reversed  and  the  same  results  achieved.   A  time- 
temperature  response  curve  of  the  fluid  stream  leaving  the  solid  is 
recorded.   By  comparing  this  record  to  theoretically  developed  non-dimen- 
sional temperature  response  curves,  the  heat  transfer  characteristics  of 
the  solid  can  be  determined. 


2.    SUMMARY  OF  THEORY. 

A  theoretical  analysis  of  the  transient  heat  transfer  behavior  of 
a  porous  medium  was  presented  by  Schumann  (2)  in  1929.   A  porous  medium 
may  be  defined  as  an  aggregation  of  solid  bodies  so  arranged  that  the 
void  spaces  between  the  bodies  are  connected  and  will  permit  the  passage 
of  a  fluid.   This  definition  can  of  course  be  extended  to  include  the 
matrix  surfaces  utilized  in  a  compact  heat  exchanger. 

In  Schumann's  analysis,  the  porous  medium  is  considered  to  be  homo- 
geneous and  initially  at  a  uniform  temperature,  t..   A  fluid  at  the  same 
temperature  is  flowing  through  the  medium.   At  a  certain  instant  the 
temperature  of  the  entering  fluid  is  assumed  to  change  to  a  higher  tem- 
perature t   .   The  problem  is  then  to  find  the  temperature  of  the  fluid 
and  the  solid  as  functions  of  time  and  position  in  the  matrix  or  porous 
media.   The  basis  of  the  analysis  are  the  theoretical  generalized  heat- 
ing curves  that  are  developed  from  an  energy  balance  on  the  matrix. 


I       A         ^s 
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1)  Properties  of  the  fluid  are  independent  of  temperature. 

2)  Fluid  flow  is  steady. 

3)  The  porous  solid  is  homogeneous. 

4)  The  thermal  conductivity  of  both  fluid  and  solid  is  infinite 
perpendicular  to  the  flow  direction. 

5)  Thermal  conductivity  of  fluid   is  zero  in  the  flow  direction. 

Heat  absorbed  by  the  solid  =   A  n$  Cs  3©  £*-* 
Heat  transferred  to  solid  by  convection  =   nb(ti  -  Cs)aX, 
Heat  transferred  from  the  fluid  by  convection  =  /V)Cp  "jT£~  0% 
Heat  transferred  from  solid  by  conduction  =   -  Ks  ns     ■>.    1  "^ 
An  energy  balance  for  the  fluid  and  solid  results  in: 

pjsc/^d*.   =kk?ph+  bUii-t^Jx  CD 

Let:   T=  generalized  time  variable  =  PP—    © 

WScs 


generalized  position  variable =  N 


tuL 


h  A 
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Making  the  above  substitutions,  equations  1)  and  2)  become: 

*£  --  (ts-t*) 

The  experimental  method  used  requires  that  the  above  assumptions 
and  the  following  boundary  and  initial  conditions  be  closely  observed. 
These  conditions  are  as  follows: 

1)  Initially  the  core  is  at  a  uniform  temperature. 

2)  At  time  t  =  0,  the  temperature  of  the  incoming  fluid  changes 
instantaneously. 


3)   System  is  adiabatic,  that  is,  only  boat  is  transferred  be- 
tween the  matrix  core  and  the  fluid, none  with  the  surroundings. 

With  the  idealizations,  initial  and  boundary  conditions,  and  the  additional 
assumption  that  the  thermal  conduction  in  the  solid  is  zero  parallel  to 
the  direction  of  fluid  flow  (i.e.,  A  ,  the  material  conduction  para- 
meter, =  0),  Schumann's  solutions  are  as  follows: 


t*  "  t. 


--(*  +T) 


=  1  -  e  v"  T  '  '      /    "    d(zT) 


,..„  j 


fi   i  / .     d(z  r  ) 


j0  (ZivTrz) 


This  solution  was  first  used  as  the  basis  for  a  transient  technique  by 

Furnas  (16)  in  1932.   The  fluid  temperature  was  measured  at  the  exit  of 

the  matrix  where  x  =  L,  therefore  z  =  N   .   The  experimental  data  was 

tu 

compared  with  Schumann's  theoretical  constant  z  curves.   The  theoretical 

z  curve  which  best  fit  the  experimental  data  was  considered  as  the  N 

tu 

of  the  matrix  and  from  this  value,  h,  the  convection  conductance  was 
determined.   The  biggest  disadvantage  of  this  method  is  in  the  calcula- 
tion of  an  infinite  series  for  each  point  and  the  number  of  terms  neces- 
sary for  a  good  approximation  increases  as  N   increases. 

In  1950,  G.  L.  Locke  (1)  developed  another  method  for  the  evaluation 

of  experimental  N   from  the  theoretical  curves.   In  this  method,  only 

tu  '     } 

the  maximum  slopes  of  the  theoretical  and  experimental  curves  need  to  be 
compared.   He  derived  an  expression  for  the  slope  of  the  generalized  heat- 
ing curve,  (t.  -  t.)/(tf.  -  t.)  vs.  Y/z,  as  a  function  of  N   .   Locke's 


expression  shows  that  the  maximum  slope  of  the  generalized  heating  curve 

is  a  unique  function  of  N   : 
tu 


d(  f2 


d(T/Ntu)     ^  tuT 


-  V2i^7  >   e"^tu+T  > 


The  above  implies  that  the  downstream  fluid  temperature  is  measured  at 

x  =  L.   Therefore,  z  =  N   and  tc  =   t-0.   This  method  has  become  known 
tu      f    f2 

as  "The  Maximum  Slope  Method". 

One  of  the  idealizations  made  by  Schumann  and  followed  by  Locke  was 
that  the  thermal  conductivity  of  the  solid  parallel  to  the  flow  is  zero* 
This  is  very  difficult  to  achieve  in  an  experiment  with  a  continuous  flow 
channel  in  the  matrix.   Although  the  effect  of  finite  thermal  conductiv- 
ity parallel  to  the  flow  is  negligible  at  high  Reynold's  numbers,  it  can- 
not be  ignored  at  low  Reynold's  numbers.   Recent  works  by  Crestwick  (3), 
Mondt  (4),  and  Howard  (5)  consider  the  effects  of  longitudinal  thermal 
conduction  in  the  solution  of  Schumann's  transient  heating  problem.   To 
correlate  results,  they  utilize  the  conduction  parameter,  A  .   This 
parameter  is  the  ratio  of  conduction  heat  transfer  to  the  change  in  in- 
ternal thermal  energy  of  the  flowing  fluid  for  the  same  temperature  differ- 
ence.  When  A  -   0,  the  maximum  slope  is  a  unique  function  of  N   s  as  put 
forth  by  Locke,  however,  when  A  f   0,  the  maximum  slope  is  a  function  of 

N   and  A  . 
tu 

Howard,  in  his  paper  (5),  included  graphs  and  a  table  of  maximum 

slopes  as  a  function  of  N   and  A   .   He  derived  these  values  by  finite 

tu  J 

difference  calculations  on  a  digital  computer  and  proposed  that  the  values 
of  maximum  slopes,  as  presented,  are  within  +  2%  for  a  given  N   and  A  . 


Ideally  then,  if  the  idealizations,  initial  and  houndary  conditions 
can  be  adhered  to,  all  that  is  required  for  the  experimental  evaluation 
of  the  heat  transfer  characteristics  of  a  matrix  is  to: 

1)  Determine  the  conduction  parameter  A  . 

2)  Experimentally  determine  the  maximum  slope  of  the  heating  or 
cooling  curve. 


3,    OBJECTIVES. 

Tn  the  past,  there  has  been  a  variety  of  experiments  conceived  to 
utilize  the  principles  previously  discussed,  by  approaching  the  idealiza- 
tions, initial  and  boundary  conditions  as  closely  as  possible.   Perhaps 
the  most  difficult  condition  to  achieve  is;  "at  time  t  =  0  the  tempera- 
ture of  the  incoming  fluid  changes  instantaneously". 

Locke  (1)  and  Coppage  (9)  both  employed  an  air  switching  device  with 
which  air  could  be  quickly  directed  through  the  matrix  or  by-passed  to 
the  atmosphere.   By  use  of  this  switch  they  imposed  a  "step  change"  in 
temperature  on  their  test  matrices.   To  approach  the  "step  change"  in 
temperature,  "requires  a  fast  operating  fluid  switch  for  introducing 
the  hotter  (or  colder)  fluid  in  order  to  assure  an  instantaneous  change 
of  fluid  temperature.   However,  even  with  rapid  switching,  one  of  the 
major  difficulties  here  is  the  inevitable  thermal  capacity  of  ductings 
switch  device  and  matrix  holding  device  causes  unavoidable  changes  in  the 
entering  temperature".  (9) 

Solar  Aircraft  Company  recently  did  some  work  testing  core  modules 
for  their  T-600  engine  (10).   Tn  their  experimental  rig,  they  utilized 
two  steady  streams  of  air  and  a  carriage  tray  that  houses  the  core  module. 
The  tray  was  inserted  through  the  side  of  the  test  cabinet  and  was  moved 
from  a  warm  air  stream  to  a  ambient  or  cold  air  stream  to  simulate  a  step 
flunction  change  of  temperature.   This  type  of  switching  was  done  manu- 
ally. 

Howard,  in  his  work  with  new  compact  surfaces  (6),  employed  a  "slid- 
ing drawer".   This  drawer  holds  the  matrix  and  the  downstream  thermo- 
couples.  The  drawer  can  be  inserted  or  withdrawn  from  the  test  section 


through  which  ambient  air  flows  continuously.   The  matrix  is  heated 
when  the  drawer  is  in  the  withdrawn  position  by  an  independent  tempera- 
ture controlled  air  stream. 

One  of  the  objectives  of  this  thesis  was  to  construct  and  test  an 
experimental  rig  that  would  improve  upon  those  previously  used.   This 
test  apparatus  described  herein  differs  from  other  test  facilities 
chiefly  in  the  method  used  to  simulate  the  initial  condition  of  a  step 
change  in  temperature  of  the  fluid  (air). 

The  test  apparatus  constructed  has  but  one  air  stream,  the  air  is 
heated  by  a  high  resistance  nichrome  wire  heater  complex,  an  integral 
part  of  the  apparatus,  located  just  upstream  of  the  test  matrix.   The 
time  constant  or  response  of  the  heaters  is  at  least  an  order  of  magni- 
tude smaller  than  the  response  of  the  thermocouples  and  recording  devices 
used.   By  turning  the  heaters  on  or  off,  the  step  change  in  temperature 
is  simulated.   The  other  systems  required  some  manual  dexterity  on  the 
part  of  the  experimenter  to  begin  a  data  run,  whereas  this  system  in- 
volves only  turning  off  or  on  the  heaters.   The  temperature  recording 
device  is  started  prior  to  the  start  of  a  run.   The  results  achieved  by 
this  system  were  an  improvement  over  other  systems.   With  the  "sliding 
drawer"  and  "quick  changing  valve"  there  is  a  change  in  the  flow  and  there- 
fore the  hydrodynamic  as  well  as  the  thermal  boundary  layer  must  develope. 
In  this  system,  the  hydrodynamic  boundary  layer  is  already  formed,  only 
the  thermal  boundary  layer  is  affected.   A  complete  description  of  equip- 
ment is  presented  in  APPENDIX  B. 

In  addition  to  the  construction  of  this  test  apparatus  and  the  verifi- 
cation of  its  reliability  by  comparing  with  published  data,  extension  of 
the  curves  published  by  Howard  (6)  into  the  high  Reynold's  number  region 
was  attempted. 


Finally,  heat  transfer  and  flow  friction  characteristics  of  several 
new  surfaces  were  determined. 
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4.    EXPERIMENTAL  TECHNIQUES 

In  reference  (5)  Howard  presents  results  from  a  numerical  finite- 
difference  method  of  calculating  the  transient  behavior  of  a  porous 
media  when  subjected  to  a  step  change  in  fluid  temperature.   These  cal- 
culations consider  the  case  where  longitudinal  thermal  heat  conduction 
cannot  be  neglected.   These  data  are  intimately  involved  in  the  determina- 
tion of  N   and  subsequently  the  Colburn  j  factor  as  presented  in  this 
thesis.   For  this  reason  the  restrictions,  idealizations  and  conditions 
Howard  placed  on  his  solution  will  be  restated  here.   It  is  these  re- 
strictions, idealizations  and  conditions  that  the  experimental  set-up 
must  approach  as  closely  as  possible  to  achieve  reliable  results. 

1)  The  flow  entering  the  matrix  is  steady  and  uniform  in  velocity 
and  temperature,  and  remains  steady  and  uniform  at  any  cross 
section  as  it  flows  through  the  matrix. 

2)  The  thermal  conductivity  of  the  matrix  is  infinite  in  the 
direction  normal  to  the  flow  and  finite  in  the  direction  para- 
llel to  the  flow. 

3)  The  thermal  capacity  of  the  fluid  contained  at  any  time  within 
the  matrix  is  small  compared  with  the  thermal  capacity  of  the 
matrix. 

4)  The  thermal  properties  of  the  fluid  and  matrix  are  constant 
and  uniform. 

5)  The  convective  heat  transfer  coefficient  is  some  suitable 
average  and  remains  uniform  and  constant. 

6)  At  time  zero  the  change  in  fluid  temperature  will  be  a  step 
function,  with  the  matrix  and  its  entrained  fluid  initially  at 
some  uniform  and  equal  temperature. 

The  first  idealization  was  attained  by  careful  experimentation. 

Uniform  velocity  profiles  were  obtained  by  wire  screen  flow  straighteners. 

Uniform  temperature  profiles  were  obtained  by  a  wire  screen  type  heater 

explained  in  detail  in  APPENDIX  B.   The  velocity  profiles  were  taken  at 

the  inlet  of  the  matrix  over  entire  range  of  flow  rates  to  ensure  that 
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the  profiles  were  uniform.   The  temperature  profiles  were  taken  over  the 
entire  range  of  flow  rates  and  found  to  be  uniform  within  +  ,5°F. 

The  second  idealization  cannot  be  strictly  realized  in  practice.   The 
effect  of  the  deviations  from  this  idealization  can  be  made  negligible  in 
the  solid  if  the  resistance  to  heat  transfer  within  the  solid  is  small 
compared  to  the  heat  transfer  resistance  between  the  solid  and  fluid. 

The  third  idealization  restricts  the  fluid  to  a  gas  rather  than  a 
liquid. 

The  fourth  idealization  can  be  realized  if  the  temperature  variation 
is  limited  to  20  -  25°F.   rise  above  ambient.   This  would  place  the 
operating  temperature  range  at  60  to  80°F.   In  this  range  the  thermal 
properties  of  the  fluid,  particularly  a  gas,  vary  insignificantly  and  can 
be  considered,  as  a  good  approximation,  constant  and  uniform. 

The  fifth  idealization  is  imposed  upon  the  computer  solution  and,  of 
course,  the  determination  of  the  convective  heat  transfer  coefficient  is 
the  ultimate  result  of  the  experimental  process. 

Finally,  two  sets  of  thermocouples  were  placed  one  upstream  and  one 
downstream  of  the  matrix.   Prior  to  the  beginning  of  each  run,  these 
thermocouple  complexes  were  read.  When  a  steady  temperature  differential 
was  indicated,  the  run  was  begun. 

Herein  lies  the  improvement  over  the  "sliding  drawer"  technique.   A 
series  of  nichrome  wire  heaters  were  installed  in  the  air  stream  to  heat 
the  air  and  subsequently  the  matrix.   To  begin  a  runs  the  heaters  were 
turned  on  or  off,  and  the  temperature  vs.  time  recorder  started.   The  time 
constant  (TC)  of  the  nichrome  wire  heater  is  very  small  (TC  =  .0425  sec 
for  m  =  1000  lb/hr,  TC  =  .259  sec  for  m  =  10  lb/hr)  and  hence  the  approach 
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to  the  final  idealization:   at  time  t  =  0  the  change  in  fluid  temperature 
will  be  a  step  function. 

A  schematic  diagram  of  the  experimental  system  in  which  air  is  used 
as  the  working  fluid  is  shown  in  Figure  1.  Figure  2  is  a  photograph  of 
the  entire  test  apparatus.   Essentially,  the  system  consists  of  a  source 
of  air,  an  air  metering  device,  flow  straightener  screen  pack,  an  air 
heater  section,  a  test  matrix,  a  pressure  measuring  system  and  a  tempera- 
ture measuring  system. 

Air  is  drawn  through  the  apparatus  by  a  six-stage  centrifugal  air  com- 
pressor.  Air  flow  rates  are  measured  by  an  ASME  standard  D  and  T)/2   ori- 
fice meter  with  changeable  orifice  plates.   Mass  flow  through  the  system 
is  controlled  by  two  sliding  valces,  one  on  the  inlet  to  the  compressor 
and  one  on  the  compressor  outlet.   In  addition,  there  is  a  fine  control 
by-pass  valve  located  just  ahead  of  the  compressor.   As  the  air  passes 
through  the  orifice  plate  it  comes  to  the  screen  flow  straightener  at 
which  point  the  flow  cross-section  changes  from  round  to  square.   The  air 
then  flows  through  the  air  heater  section  which  serves  a  dual  purpose  of 
straightening  the  air  velocity  profile  and  heating  the  air.   In  this 
section  the  air  is  heated  uniformly  to  a  predetermined  temperature,  usually 
20°F  above  ambient.   The  heated  air  passes  into  the  matrix  and  heats  the 
matrix  to  the  same  uniform  temperature.   At  this  point,  a  run  is  ready  to 
start.   By  referencing  the  thermocouples  downstream  of  the  matrix  to  the 
thermocouples  upstream  of  the  heaters,  the  initial  temperature  differen- 
tial can  be  very  closely  controlled.   This  temperature  differential  is 
continuously  recorded  by  a  Leeds  and  Northrup  "Speedomax"  strip  chart 
recorder  when  a  run  is  started.   The  recorded  trace  has  a  distinct  ad- 
vantage in  that  there  is  no  transposition  of  data  necessary,  hence,  a 
decrease  in  the  uncertainty  of  the  results.   The  resultant  trace  or  curve 
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is  the  primary  data  for  determining  j.   The  N   of  the  matrix  can  be  evalu- 
ated by  determining  the  maximum  slope  of  this  fluid  temperature  differ- 
ence curve  during  the  cooling  transient.   This  is  the  basis  of  the  maxi- 
mum slope  technique  as  developed  by  Locke  (1), 

The  flow  friction  performances  of  the  matrices  tested  were  obtained 
by  determining  the  isothermal  pressure  drop  across  the  matrix  and  the 
mass  flow  rate  of  the  fluid  through  the  matrix.   The  pressure  drop  data 
were  obtained  from  static  pressure  taps  located  in  the  test  section  im-  , 
mediately  upstream  and  downstream  of  the  test  section.   The  entering 
fluid  static  pressure  was  also  recorded.   The  mass  flow  rate  data  were 
obtained  in  the  air  metering  section  as  previously  described.   In  summarys 
the  following  data  were  recorded  for  each  run: 

a)  Orifice  diameter  (D  ) 

b)  Pressure  drop  across  the  orifice  (  AP  ) 

c)  Pressure  drop  across  the  matrix  (APT 

d)  Entering  static  pressure  (Ps) 

e)  Orifice  temperature  (T  ) 

f)  Atmosphere  pressure  (P  ) 

g)  Temperature-Time  cooling  or  heating  curve 
h)  Temperature  recorder  chart  speed 

A  complete  description  of  the  experimental  apparatus  and  equipment  is 

in  APPENDIX  B.   APPENDIX  A  discusses  the  relation  of  the  data  to  the 

equations  governing  j  &  f. 
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5.    DESCRIPTION  OF  SURFACES. 

There  are  data  presented  in  this  thesis  for  seven  surfaces.   Of 
the  seven,  three  are  new.   Three  of  the  remaining  four  surfaces  were 
tested  by  Howard  (6).   The  remaining  surfac   was  used  primarily  as  a 
test  surface  to  check  the  reliability  of  the  experimental  apparatus. 

The  test  surface  was  a  matrix  consisting  of  20  randomly  stacked 
stainless  steel,  60  x  60  mesh  creens.   A  photograph  of  the  matrix 
along  with  a  tabulation  of  its  physical  and  geometrical  properties  may 
be  found  in  Figure  3.   The  heat  transfer  and  flow  friction  characteris- 
tics for  such  a  matrix  were  published  by  Coppage  in  1952   9),  and  was 
considered  a  good  check  on  the  performance  of  the  experimental  set-up. 

The  three  matrices  examined  by  Howard  and  re-examined  in  this  thesis 
were: 

1)  Plate-fin  reference  matrix 

2)  10°  skew  passage  matrix 

3)  20°  skew  passage  matrix 

These  particular  surfaces  were  chosen  to  provide  an  additional  check 
on  the  reliability  of  the  experimental  setup  and  to  provide  a  compari- 
son between  the  "integral  heater"  and  "sliding  drawer"  experimental  techni- 
ques. 

The  plate-fin  surface  is  formed  by  corrugated  sheets  separated  by  a 
flat  spacing  sheet  to  prevent  meshing  of  the  corrugations.  This  matrix 
was  to  serve  as  a  reference,  since  the  theoretical  laminar  flow  solutions 
for  both  heat  transfer  and  frictional  behavior  of  triangular  cross-sec- 
tions are  available.  However,  difficulty  was  encountered  in  forming  the 
material  into  triangular  cross-sections,  and  therefore  it  fails  as  a  re- 
ference, but  the  results  were  interesting  and  will  be  discussed  later. 


Figure  4  presents  photographs  of  the  plate-fin  matrix  and  a  tabulation 
of  its1  physical  and  geometrical  properties. 

The  skewed  passage  also  employs  corrugated  sheets,  but  unlike  the 
plate-fin,  the  spacer  plate  is  not  included.   To  prevent  meshing  of  the 
corrugations,  the  sheets  are  skewed  relative  to  one  another  forming  what 
is  defined  as  the  "skew  angle".   A  photograph  of  the  matrix  and  a  tabula- 
tion of  its  properties  may  be  found  in  Figure  5. 

Of  the  three  new  surfaces  tested,  one  is  a  version  of  the  previous- 
ly tested  20°  skew  matrix  using  perforated  or  slotted  nickel  plates  for 
a  material  vice  a  continuous  material.   The  plates  were  corrugated  on 
the  same  device  used  for  the  orig-nal  20°  skew  to  maintain  the  same  corru- 
gation geometry.   The  perforations  should  reduce  stagnant  areas  in  the 
matrix,  thereby  increasing  the  heat  transfer  characteristics. 

Finally,  two  matrices  of  a  parallel  plate  design  x^ere  constructed 
and  tested.   The  "parallel  plate"  matrix  is  similar  to  the  plate-fin 
in  that  a  spacer  plate  is  placed  between  the  formed  sheets  to  prevent  mesh- 
ing.  The  sheets  are  formed  by  rolling.   When  assembled,  the  matrix  is  a 
maze  of  rectangular  passages  with  an  aspect  ratio  of  77.   Aspect  ratio  is 
defined  in  reference  (8)  as  the  ratio  of  the  two  sides  of  a  rectnagular 
flow  passage,  b/a.   The  sketch  below  is,  ideally,  what  was  intended  in 
designing  the  "parallel  plate"  matrix. 


One  matrix  was  formed  from  a  brass  material  and  is  used  as  a  refer- 
ence to  theoretical  laminar  flow  data  for  parallel  configurations  (8). 
The  second  matrix  was  formed  from  the  perforated  nickel  plate  to  examine 
the  effects  of  the  perforations  on  this  type  of  geometry.   Pictures  and  a 
tabulation  of  properties  for  the  "parallel  plate"  matrices  are  in  Figure  6. 
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6.    PRESENTATION  OF  RESULTS. 

The  heat  transfer  and  flow  friction  data  for  each  matrix  are  pre- 
sented in  Tables  I  to  VII.   These  data  are  also  presented  in  graphic 

2/3 

R   ; 

factor,  f,  are  plotted  as  functions  of  Reynold's  Number,  N  ,  in  Figures  7 
through  16.   Plotted  also  in  Figures  13  through  15,  for  comparison  pur- 
poses, are  the  theoretical  laminar  flow  solutions  for  parallel  plates. 
For  the  heat  transfer  data,  Prandtl  number  (N  )  was  not  really  a 

test  variable,  since  the  working  fluid  was  air  at  moderate  temperatures. 

2/3 
N      is  included  as  an  approximation  to  the  effect  of  Prandtl  number 

over  a  moderate  range  of  Prandtl  numbers.   Basic  laminar  boundary  layer 

solutions  indicate  that  in  the  Prandtl  number  range  of  0.5  to  1.0  (gases), 

the  Prandtl  number  enters  the  solution  as  approximately  the  2/3  power. 

Since  a  moderate  temperature  range  is  used  in  the  experiment,  air 
properties  throughout  the  matrix  are  practically  constant.   Therefore, 
all  the  properties  of  air  used  in  the  calculations  of  the  heat  transfer 
and  friction  characteristics  were  those  of  air  at  the  bulk  temperature  of 
the  air  stream  as  measured  at  the  orifice  inlet. 

The  Reynold's  number  is  based  on  the  hydraulic  diameter  and  mass  flov; 

velocity. 

DH  G 
NR  ~~   /X 
The  hydraulic  diameter  is  defined  as: 


n     a    -  L(     Fl°w  Passage  Cross-sectional  area 
H      H       Flow  Passage  Wetted  Perimeter     ' 

-   4   Void  Volume   x  L 

L  Heat  Transfer  Area 


The  mass  flow  velocity  is  based  on  the  matrix  free  flow  ctoss-sectional 


A  pA. 

c  v    fr 

The  friction  factor  includes  the  entrance,  exit  and  flow  accelera- 
tion effects.   The  entrance  and  exit  effects  normally  provide  only  a 
small  contribution  to  the  overall  pressure  drop  in  a  matrix.   Core  fric- 
tion is  the  term  that  controls  the  magnitude  of  the  pressure  drop.   As  a 
result,  high  accuracy  in  the  evaluation  of  the  entrance  and  exit  loss  co- 
efficients is  not  required.   These  coefficients  k  and  k  are  functions 

c      e 

of  flow  geometry  and  Reynold's  number  and  were  obtained  from  Chapter  IV, 
Figures  18  to  21,  as  appropriate,  reference  (8). 
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7.    UNCERTAINTIES  AND  LIMITATIONS. 

As  put  forth  in  the  theory,  the  experimental  method  utilized,  re- 
quires great  care  to  approach  the  idealizations  and  initial  and  boundary 
conditions  of  the  theoretical  analysis.   The  magnitude  of  errors  intro- 
duced if  the  experimental  system  does  not  duplicate  the  ideal  system  is 
difficult  to  evaluate.   Each  of  the  conditions  and  idealizations  were 
carefully  considered  and  the  experimental  system  was  designed  to  approach 
as  closely  as  possible  the  ideal  system.   These  items  are  discussed  in 
the  "Experimental  Technique"  section.   No  attempt  will  be  made  here  to 
evaluate  or  estimate  this  type  of  error. 

Other  experimental  errors  are  caused  by:  1)  uncertainty  in  the 
value  of  physical  constants,  2)  inaccuracy  in  geometric  measurements 
and  3)  instrument  inaccuracies.   The  method  employed  herein  to  describe 
the  uncertainty  in  the  values  reported  is  explained  fully  in  Reference  (14). 
However,  before  stating  these  uncertainties,  an  example  of  the  process 
follows: 

A  theorem  of  statistics  as  stated  in  Reference  (14)  is: 

"If  R  is  a  linear  function  of  n  dependent  variables,  each  of  which 
is  normally  distributed,  then  the  relation  between  the  interval  for  the 
variables,  w^,  and  the  interval  for  the  result,  wR,  which  gives  the  same 
odds  for  each  of  the  variables  and  for  the  result  is: 

-ife 
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where   w.  =  uncertainty  interval  based  on  certain  odds. 
Consider  the  Colburn  j  factor 


i  .  .       *e  -  2'3 


by  the  theorem  and  letting  P  =  N, 


tu  A   PR 
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if  both  sides  are  divided  by  j  then: 


(1) 


Before  an  equation  such  as  this  can  be  evaluated  the  uncertainty 
terms  must  be  known. 

The  physical  constants  necessary  are  listed  below.   References  (12) 
and  (13)  were  sources  for  the  values  of  these  constants.   The  uncertainties 
in  these  values,  as  closely  as  can  be  determined,  are  also  listed  below, 
c  +  0.5% 


P 


±  0.5% 
+  1.0% 
+  2.0% 
+  0.5% 


The  uncertainty  in  physical  linear  measurements  could  be  consider- 
ed negligible  if  core  construction  did  not  have  minor  inconsistencies. 
The  uncertainty  in  the  weight  of  a  core  is  considered  negligible.   Conser- 
vative estimates  of  the  uncertainty  in  linear  measurements  are; 


fr' 


L   +  0.5% 
W   <  0.1% 
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Instrumentation  errors  were  such  that  only  the  manometers  need  be 
considered.   With  the  exception  of  the  orifice  temperature,  all  tempera- 
tures measured  were  temperature  differences  and  these  were  measured  in 
inches  as  indicated  in  Appendix  A.   Since  pressures  varied  and  manometers 
to  measure  them  changed,  the  maximum  uncertainty  encountered  was  used  in 
the  error  analysis  and  is  listed  below. 

P         +  1.257. 
o        — 

P        +  1.707. 

m        — 

P  negligible 

atm 

Another  uncertainty  which  requires  special  mention  is  that  in  N 
Since  j  does  not  vary  linearly  with  N  ,  the  uncertainty  at  high  and  low 
N  was  considered.   This  uncertainty  was  determined  by  utilizing  Figure  19 
and  the  uncertainty  in  Maximum  Slope  as  2.07.. 

N   =  3.5,  X=  0 


N 
N 


30,   X=  .04   dNtu  =  +  9.67. 


Returning   to  Equation    (1)    and   substituting  as   appropriate  : 
For  high  NR:       W^  s    Loi))^  (ooV(.  0lf+  (^f  1   **■    •  ofc3 

r      *  *n£ 

For  low  NR:   w^  -      (.0%)  4  (.01)  +(.  01)  +  (;£i)    -  •  °1& 


0 


J 


Following  the  above  type  of  analysis,  the  following  additional 
uncertainties  were  calculated. 

m  +1.07. 
N  +  2.57. 
f    +  4.37. 
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8    DISCUSSION  OF  RESULTS. 

The  60  x  60  mesh  screen  matrix  was  used  solely  as  a  test  matrix,  a 
matrix  whose  heat  transfer  and  friction  characteristics  were  well  docu- 
mented.  Figure  7  and  Table  I  present  the  results.   There  is  quite  good 
agreement  in  the  friction  factor  but  a  slight  variation  in  j  factor.   This 
difference  is  attributed  to  the  difference  in  the  geometries  of  the  test 
matrices.   Exact  reproduction  of  a  screen  matrix  is  difficult  to  achieve. 
The  correlation  was  such  that  the  experimental  test  setup  was  considered 
capable  of  delivering  reliable  data. 

For  further  verification  on  the  reliability  of  the  test  setup,  Figures 
8  to  10  and  Tables  II  to  IV  present  the  test  results  for  the  Plate-Fin 
Reference  Matrix,  the  10°  Skew  Stainless  Steel  Matrix  and  the  20°  Skew 
Stainless  Steel  Matrix  used  by  Howard  (6).   Although  this  was  not  the 
original  purpose  for  utilizing  these  matrices,  the  results  do  indicate 
conclusively  that  the  "integral  heater"  is  as  good  or  better  than  the 
"sliding  drawer"  technique. 

The  original  purpose  for  utilizing  these  matrices  was  to  try  and  ex- 
tend Howard's  curves  into  the  higher  Reynold's  number  region  (lox^er  N 
region).   The  points  obtained  in  this  region  are  plotted  on  the  curves  and 
the  perinent  data  are  recorded  in  the  Tables.   None  of  them  follow  the 
straight  line  of  their  predecessors.   There  are  at  least  two  reasons  for 
this  obvious  departure  from  the  norm.   First,  it  is  a  well  known  fact 
that  the  "Maximum  Slope"  Method  has  an  inherently  large  error  in  the  low 
N   region  of  1.5<N  <3.5.   This  fact  was  originally  published  by  Locke 
(1)  for  A=  0.   Howard  (5)  took  into  account  longitudinal  conduction  and 
published  an  error  curve  that  includes  the  effect  of  longitudinal  conduc- 
tion on  an  error  in  N   ,  see  Figure  19.   In  the  regions  below  N   <   1.5 
tu'       °  6  tu 
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the  error  is  acceptable,  yet,  a  discrepancy  still  exists.   This  points 
out  the  second  reason  for  the  discrepancy.   The  maximum  slope  at  high 
flow  rates  (low  N  )  occurs  at  time  t  =  0+»   Instrumentation  prevents 
accurate  measurement  of  the  true  maximum  slope,  at  this  time  due  to  in- 
herent delay  or  lag  time. 

Extrapolation  of  the  time-temperature  curves  back  to  time  t  =  0  was 
attempted  in  an  effort  to  remove  the  effects  of  instrumentation.   The 
method  used  was  either  incorrect  or  unreliable.   Satisfactory  results 
were  unobtainable. 

The  merits  of  the  plate-fin,  10°  skew  and  20°  skew  matrices  as  heat 
transfer  surfaces  was  discussed  thoroughly  by  Howard  in  reference  (6). 

Figures  11,  13  and  14  and  Tables  V,  VI  and  VII  present  the  results 
for  a  new  geometry  tested  and  a  new  idea  in  material.   Figure  11  and 
Table  V  present  the  results  for  a  20°  skew  nickel  matrix  constructed  from 
slotted  plates.   Appendix  C  explains  fully  the  geometry  of  the  slots  and 
the  calculations  that  determine  the  basic  geometric  parameters  p,  p    , 
and  r  .   The  effects  of  longitudinal  conduction  was  considered  negligible 
as  the  matrix  was  constructed  so  that  the  air  flow  was  normal  to  the  long 
dimension  of  the  slots.   Plotting  all  the  points  uncorrected  for  conduc- 
tion showed  an  erratic  departure  from  the  expected  results  for  N   <  200. 
It  was  then  decided  that  perhaps  the  basic  assumption  was  not  warranted 
completely.   Corrections  for  longitudinal  conduction  were  made  where 
significant  and  these  points  were  plotted.   These  corrected  values  of  j 
are  overestimations,  for,  although  A  is  corrected  for  the  effects  of  the 
slots,  the  matrix  length  L  is  not  the  true  conduction  path  length  but  at 
best  a  lesser  approximation.   Since  A  is  obviously  not  equal  to  zeros 
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the  true  value  of  j  is  probably  between  the  uncorrected  j  (  A=  0)  and 
the  j  corrected  for  longitudinal  conduction  using  L  as  an  approximation 
for  conduction  path  length.   This  was  considered  in  the  interpretation  of 
results. 

Figure  12  compares  the  20°  skew  stainless  steel  matrix  with  the  per- 
forated nickel  matrix.   There  is  a  70  to  1007o  increase  in  j,  while  f  in- 
creases 10  to  25%. 

It  must  be  understood  that  these  two  matrices  are  geometrically  similar 
only  in  overall  dimensions.   Actually,  the  heat  transfer  areas  were  cal- 
culated on  different  bases.   The  20°  skew  stainless  steel  heat  transfer 
area  was  based  on  the  total  surface  area  available  for  heat  transfer  while 
the  heat  transfer  area  of  the  20°  skew  perforated  was  based  on  this  same 
area  less  the  area  of  the  slots.   Also,  not  included  in  the  heat  trans- 
fer area  is  the  peripheral  area  of  the  slots  normal  to  the  flow.   If 
the  20°  skew  perforated  j  factor  were  based  on  the  total  available  area 
for  heat  transfer,  the  values  plotted  would  be  decreased  approximately 
40%.   The  perforated  material  still  represents  a  substantial  improvement 
in  the  performance.   This  improvement  is  attributed  to  the  reduction  of 
inefficient  stagnant  areas  within  the  matrix  and  the  alteration  of  the 
thermal  and  hydrodynamic  boundary  layers  caused  by  the  slots. 

Figure  13  and  Table  VI  present  the  results  for  a  "parallel  plate"  mat- 
rix constructed  from  brass  shim  stock.   As  can  be  seen  in  Figure  13,  the 
heat  transfer  performance  falls  far  below  that  of  theoretical  laminar 
flow  performance  for  parallel  plates.   Several  things  may  account  for  this: 
the  aspect  ratio  of  the  flow  channels  in  this  matrix  was  77  vice  CO     for 
truly  parallel  plates,  and  also,  the  flow  channels  were  not  uniform  in 
cross-section  from  one  to  another.   The  non-uniformity  in  cross-section 
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results  in  uneven  mass  flow  through  the  matrix,  hence  inefficient  use 
of  the  available  heat  transfer  surfaces,  and  therefore  a  decrease  in  j. 

The  decrease  in  friction  factor  is  possibly  due  to  the  sharp  angle 
regions  where  the  formed  sections  contact  the  flat  plate.   The  wall  shear 
stress  in  this  type  of  region  is  small  which  would  tend  to  reduce  the 
value  of  f.   These  regions  would  also  tend  to  produce  stagnant  regions 
previously  mentioned,  having  a  diminishing  effect  also  on  j„ 

Figure  14  and  Table  VII  present  the  results  for  a  "parallel  plate" 
matrix  constructed  from  perforated  nickel  plate.   This  matrix  was  construct- 
ed in  the  identical  fashion  as  the  brass  matrix.   Although  it  is  not  read- 
ily apparent  in  Figure  14,  the  heat  transfer  characteristics  have  been 
improved  nearly  807o,  while  the  friction  factor  has  remained  the  same. 
Again  it  is  necessary  to  make  note  of  the  fact  that  the  heat  transfer  area 
of  the  two  matrices  were  calculated  on  different  bases.   These  bases  are 
the  same  as  discussed  for  the  20°  skew  matrices.   Applicable  also  is  the 
approximate  reduction  by  40%  of  the  plotted  values  of  the  j  factor  if 
all  of  the  area,  in  the  perforated  material,  available  for  heat  transfer,, 
were  considered  in  the  calculations. 

Figure  15  is  a  comparison  of  the  two  "parallel  place"  matrices.   The 
improvement  realized  by  using  perforated  material  is  readily  apparent, 
even  considering  the  40%  reduction  mentioned  above.   The  additional  per- 
formance gain  in  the  transferral  of  heat  is  again  attributed  to  the  re- 
lieving of  the  stagnant  areas  and  the  alteration  of  the  boundary  layers 
caused  by  the  slots. 

The  fact  that  the  friction  factor  did  not  increase  does  not  lend 
itself  to  previous  conclusions  (20°  skew);  that  the  presence  of  slots 
increases  the  friction  factoi .   Some  additional  factor  must  be  necessary. 
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tps,  the  skewing  coupled  with  the  slots  contribute  more  to  friction 
than  straight,  slotted  channels. 

Tn  the  determination  of  the  j  factor  for  the  perforated  "parallel 
plate"  matrix  the  conduction  parameter,  A  ,  was  considered  zero.   Here 
again  the  plotted  points  departed  from  the  expected  straight  line  indicat- 
ing that   A  >  in  fact,  was  not  zero.   Corrections  for  longitudinal  conduc- 
tion were  made  using  the  matrix  length,  L,  as  an  approximation  for  the 
conduction  path  length.   These  points  were  plotted  also.   Since  the  cor- 
rected points  are  overestimat  ions  and  a/  0,  the  true  value  of  j  lies 
between  these  limits. 

The  designer  of  a  compact  heat  exchanger  is  generally  restricted  as  to 
frontal  area,  volume,  weight  and  pressure  drop.   There  is  presently  coming 
into  use,  several  means  of  comparing  heat  transfer  surfaces  with  these 
restrictions  in  mind.   Essentially,  various  parameters  are  compared 
graphically  so  that  the  relative  merit  of  the  various  surfaces  may  be  es- 
tablished. 

The  ratio  j/f  is  plotted  vs.  Reynold's  Number  as  a  means  of  relating 
frontal  flow  area  with  pressure  drop.   For  a  given  pressure  drop  a  high 
j/f  characteristic  implies  a  small  flow  frontal  area.   j/f  when  plotted 
vs.  N   is  known  as  the  Flow  Frontal  Area  "Goodness"  factor. 

The  heat  transfer  power,  h    ,  is  plotted  vs.  the  flow  friction  power, 
E    ,  and  the  resultant  characteristic  curve  is  a  measure  of  relative  core 


known  as  the  Volume  "Goodness"  factor. 

Figure  16  compares  the  20°  skew  matrices  and  the  "parallel  plate" 
matrices.   Although  the  "Goodness"  factors  are  not  infallible,  they  are 
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indicative  of  relative  performance.   In  this  light,  the  perforated 
material  should  be  excellent  matrix  material. 
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9.    CONCLUSIONS. 

Based  upon  the  correlation  of  present  data  with  that  obtained  by 
previous  experimenters,  the  experimental  test  rig  design  presented  here- 
in is  considered  capable  of  producing  reliable  test  data  on  the  heat  trans- 
fer and  friction  characteristics  of  a  core  module. 

Basic  heat  transfer  and  flow  friction  characteristics  for  a  new 
geometry  and  a  new  idea  in  material  are  presented  in  Figures  11  to  15. 
The  geometric  and  physical  properties  are  summarized  in  Figures  5  and  6. 
There  appears  to  be  a  definite  advantage  to  constructing  matrices  from 
perforated  material.   Performance  is  improved  considerably  and  based  on 
the  "Goodness"  factors  (Figure  16)  size  and  weight  advantages  are  realiz- 
ed. 

The  20°  skew  perforated  nickel  matrix  was  the  top  surface  in  perfor- 
mance.  This  matrix  had  the  highest  Flow  Frontal  Area  "Goodness"  factor, 
as  well  as  the  highest  Volume  to  surface  "Gcodness"  factor.   This  does  not 
mean  that  this  is  the  "optimum"  surface  for  any  given  application  but  it 
does  merit  consideration. 

In  both  geometries  tested  (20°  skew  and  "parallel  plate")  the  per- 
forated material  out-performed  its  counterpart.   This  would  also  indicate 
that  further  testing  of  this  type  of  material  (slotted)  is  warranted. 
It  would  seem  that  this  material  and  some  geometry  should  result  in  an  out- 
standing heat  exchanger  surface. 


10,   RECOMMENDATIONS  FOR  FUTURE  WORK. 

1)  Some  means  should  be  found  to  reliably  determine  N   in  the 

tu 

high  N   range  (1.0  ^  N   <"  3.5).   Most  promising  is  the  use 

of  a  periodic  or  cycling  technique  in  which  effectiveness  is 

measured  and  used  to  find  N   rather  than  the  maximum  slope, 
tu  t 

This  ability  is  necessary  if  performance  curves  are  to  be  extend- 
ed to  N  *>    500. 

2)  More  work  can  be  done  with  perforated  material.   Different 
geometries  should  be  tested  to  see  if  consistantly  high  per- 
formances are  realized  as  they  have  been  in  the  past.   Addition- 
ally, this  type  of  surface  is  suspect  of  fouling.   If  these 
materials  are  as  good  as  they  seem  to  be,  then  examination  of 
the  effects  of  fouling  on  their  performance  is  in  order. 

3)  As  more  surfaces  become  available,  they  should  be  tested 

and  their  heat  transfer  and  friction  characteristics  determined. 
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18-8  St. 
Steel 


Matrix  Material 

Specific  Heat  (cs)  BTU/1d  °F  °-l2 

Thermal  Conductivity  (kg)  BTU/hr  •*  ft  10.0 


Material  Thickness,  inches 
Total  Heat  Transfer  Area  (A)  ft 
Frontal  Area  (Afp)  ft 
Total  Conduction  Area  (Ag)  ft 
Plow  Cross  Section  Area  (Ac)  ft 
Matrix  Volume  (Vm)   ft 
Matrix  Density  (yOm)   lb/ft 

Hydraulic  Diameter  (3>H)  ft 
Compactness  (y9)  ft  /ft 
Porosity  (p) 


0.017 

4.0159 

0.07055 

0.02215 

0.0484 

0.00199 

155. 

0.00136 

2009. 

0.686 


FIGURE  3  Physical  and  Geometric  Properties  of 
60x60  Mesh  Screen  Matrix 
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Matrix  Material 

Specific  Heat  (cg)  BTU/lb  °F 

Thermal  Conductivity  (k  )  BTU/hr  °F  ft 


Type  302 
S.  Steel 

0.11 

7.0 


Material  Thickness,  inches 

Total  Heat  Transfer  Area  (A)   ft2 

Frontal  Area  (Afr)  ft2 

Total  Conduction  Area  (A0)  ft 

2 
Flow  Cross  Section  Area  (A  )  ft 

Matrix  Volume  (V  )  ft5 

Matrix  Density  fo)   lb/ft5 


Hydraulic  Diameter  (D„) 
Compactness  (/9)  ft2/ft5 
Porosity  (p) 


ft 


0.001 
30.12 
0.07023 
0.00752 

0.06271 
0.01170 
52.9 

.001388 

2573. 

0.893 


FIGURE  4 


Physical  and  Geometric  Properties  of 
Plate-Fin  Reference  Matrix 


SKEW   20°  SKEW 


Matrix  Material 

Specific  Heat  (V)  BTU/lb  °F 

Thermal  Cond.  (kg)  BTTJ/hr  °F  ft 

Material  Thickness,  inches 
Total  Heat  Transfer  Area  (A)  ft 
Frontal  Area  (Afp)  ft2 

Total  Conduction  Area  (Ae)  ft 

2 
Flow  Cross  Section  Area  (A  )  ft 

Matrix  Volume  (Vm)  ft3 

Matrix  Density  (0)  lb/ft5 

Hydraulic  Diameter  (DR)  ft 
Compactness  (/3)  ft  /ft5 
Porosity  (p) 


FIGURE  5  Physical  and  Geometric  Properties  of 
Skew  Matrices 


Type  302 
S. Steel 

Type  302 
S. Steel 

perf. 
Nickel 

0.11- 

0.11 

0.1065 

7.0 

7.0 

38.7 

0.001 

0.001 

0.0022 

18.38 

18.04 

15.41 

0.07024 

0.0702 

0.0734 

0.00464 

0.00452 

0.00349 

0.0656 

0.0657 

0.06412 

0.0117 

0.01157 

0.01211 

32.4 

32.22 

50.1 

.002381 

.002416 

.002528 

1570. 

1549. 

1272. 

0.934 

0.936 

0.873 

■, 

mi  i|i'|iiri|i|i|ifi|i|in 

Matrix  Material 

Specific  Heat  (cg)  BTTJ/lh  °F 

Thermal  Conductivity  (>  )  BTU/hr  °F  ft 


Material  Thickness,  inches 

Total  Heat  Transfer  Area  (A)  ft 

Frontal  Area  (A~  )  ft 

? 
Total  Conduction  Area  (A„)  ft 

2 
Flow  Cross  Section  Area  (A  )  ft 

Matrix  Volume  (Vm)  ft5 

Matrix  Density  S)   lb/ft5 

Hydraulic  Diameter  (D„)  ft 
Compactness  (0)  ft  /ft^ 
Porosity  (p) 


70-50 
Brass 

Perf. 
Nickel 

0,092 

0.1065 

66. 

38.7 

0.001 

0.0022 

22.057 

18.586 

0.06986 

0.0695 

0.00551 

0.00416 

0.06434 

0.0584 

0.01164 

0.01158 

41.5 

63.6 

.001944 

.001928 

1895. 

1605. 

.921 

.840 

FIGURE  6 


Physical  and  Geometric  Properties  of 
"Parallel  Plate"  Matrices 
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Figure  7  60x60  Mesh  Stainless  Steel  Screens 


Figure  8  *  Plate  Pin  Reference  Matrix 


1.0 


.001 


Figure  9  10°  Skew  Stainless  Steel  Matrix 


1.0 


Figure  10  20°  Skew  Stainless  Steel  Matrix 


Figure  11  20°  Skew  Nickel  (Perforated)  Matrix 


1.0 


Yftf.'m*-&*- 


Figure  12  Comparison  of  perforated 

and  non-perforated  20°  Skew 

.  Matrices 


Figure  13     "Parallel  Plate"  Brass  Llatrix 


Figure  14  "Parallel  Plate"  Nickel  (perforated)  Matrix 


1.0 


Figure  15  Comparison  of  Perforated  and 

non-perforated  "Parallel  Plate" 
Matrices 
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(2) "Parallel  Plate"  Perforated  Nickel 
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0  "Parallel  Plate"  Brass 


Figure  16  Design  "Goodness"  Factors 


TIME 

FIGURE   17    GENERALIZED    COOLING     CURVE 

NOTE :  All  measurements   are    made    in  inches 
to    facilitate    computations. 
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FIGURE  22.  Matrix  Holder  modified  with  probe 
for  velocity  profile  measurement 
and  temperature  distribution 
measurement  (Downstream  View) 


FIGURE  23.  Nichrome  Wire  Heater  (Upstream  View) 


FIGURE  2Uo  Nichrome  Wire  Heater  (Downstream  View) 


r-i  CV!  C^-CnAsO   C-OO    OS 


1 

[;■; 

' 

\ 

i, 

_L_  ./ 

■- 

-1'  !       !    ;      < 

> 

\  l 

A 

.!l|\ 

I'ii- 

f 

!( 

;|| 

to 

p 

H 

i 

ml 

| 

|k 

^ 

i  -■!' 

\    l'i« 

ill 

,  ilwil  /   i  ii 

\ 

!;: 

II 

i!  ■ 

j 

iW  ^J  -'K\            HH  i 

1 

! 

I 

i  i 

' 

:     1     : 

■  i|f 

, 

\H 

\i 

/ 

1 

:      ■   ' 

- 

1Mb 

j  \s 

r 

/ 

/ 

J\ 

IIP 

/ 

y 

\  T 

HUH 

i ! 

/ 

/ 

Xs 

\ 

•  i» 

\i 

ilit 

\ 

K>i'hllH         tw 

\\  i 

i 

III 

hi/ 

i||  I  T 

1 

jliif  = 

iii/ 

llll 

,/ 

illl 

3 '  \.  I             6 ' 

1 

J"' 

■ 

3 

\ 

"v         N 

jj;"4| 

l!!  1 

•■1  ' 

m  ii 

lljl 

■ii: 

\ 

\  l! 

/ 

1     1 

iiii|nii 

ill 

HI 

;:; 

N 

\ii 

Nil 

Ml 

1 

■  1 

! 

' 

■ 

\ 

,  V 

i 

1 

ii 

l!!i 

Ii 

\ 

|lh 

\l' 

:N 

1 

i  i1 

/ 

: 

III 

|UJ 

V 

: 

III 

"\ 

| 

% 

U|l 

1 

m 

? 

rp 

oI-j|      ^       „. 

\ 

\\\ 

Ms 

\' 

! 

< 

,. 

:! 

H 

* 

ai         ^ 

P     \ 

\ 

\ 

III 

------ - 

hi 

i 

jjfjjgjr 

|i 

\ 

Vi 

I 

=  r 

{1*1  i 

■ 

l 

I'll 

|  i     Sj        a 

1      PI* 

III! 

1 

i: 

I!,. 

I 

II  liiL 

1 

1    |i||,liiii| 

6                         '71 ;  j ' 1 

3 

i  ||l 

ii 

3'  !  '             to 

Ill' 

|| 

mi 

111 

1 

lli!l 

i     :i! 

iHliiiiblilHl 

fi'l 

!:! 

!jj! 

II       ! 

I!  ||||| 

nil1!  |H 

' 

1 

'■•illli 

; 

! 

-     i 

TT 


PROGRAM  MAXSLO 
L  =  2 
K=1 
30  READ  300,XML.AS,RH,POR,AFL,AHT,  WM 
300  F0RMAT(7F10.5) 

40  READ  400, EK,CK,SK,CM,N 
400  F0RMAT(4F10.5,  12) 

m=i  .:■•;-. 

10  READ  100,D0,DELP0,GAM, ATMP,COR, UH, UFR 

100  FORMAT(7F10.5) 

READ  101,CS,DELH  , HS , SLO, TEMPO 

101  F0RMAK5F10.5) 
EMDOT=359.0»DO*DO«SQRTF(DELPO»GAM*ATMP/29.92)*COR 
CP=0.2U 

RNUMP=5.09»EMD0T/UFR 

RNUMH=(4.0»XML*EMD0T)/( AHT*UH)  ./'• 

RNUMFR=RNUMH«UH/UFR 
CAPFL=EMDOT»CP/ 3600.0 
CAPS=WM*CM 

SLOMAX=(CAPS»SLO)/(CAPFL»CS) 
CONDPAR=(SK*AS)/(XML*EMCOT*CP) 
B=CK+EK 
C=1.0+POR«*2 
G=EMD0T/(AFL»3600.0) 
DELP=DELH». 03613 
PS=HS«. 03613 
PA=ATMP».4912 
P1=PA-PS 
P2=P1-DELP 
PM=(Pl+P2)/2. 

RHOM=(PM»144.0)/(53.3»(TEMPO+460.0)  ) 
FF1=( (64.4»RH0M»DELP*144.0)/(G»C) ) 
FF2=C»DELP/PM 
FF3=FFl-FF2-8 
FFR=FF3»RH/XML 
WRITE    OUTPUT    TAPE    4,4, M 
4    F0RMAT(9H    RUN    NO.     12///) 

WRITE    OUTPUT    TAPE    4, 1 , EMDOT , RNUMP.RNUMH, CAPFL.CAPS 

1  F0RMAT(7H    EMD0T=F1 0.5 , 3X ,7H    RNU^P=E 20.8 , 3X , 7H    RNUMH*F10 . 5, 3X, 
17H    CAPFL=F10.5,3X,6H    CAPS=F1 0.5///) 

WRITE    OUTPUT    TAPE    4, 2,CGNDPAR, S LOMA X 

2  F0RMAT(9H    C0NDPAR=F1 0.5 , 3X ,8H    SLOMAX=F 1 0. 5/// ) 
WRITE    OUTPUT    TAPE    4, 3, G , RNUMFR, FFR 

3  F0RMAT(3H    G=F10.5, 3X, 8H    RNUMFR=F10.  5, 3X, 5H    FFR=F10.5/> 
M=M*1 

IF(N-M)50, 10,10 
50    K*K+1 

IF(L-K)20f30,30 
20    STOP 

END 

END  _  „ 


Figure   27     Program  MaxSlo   Readout- 


TABLE  I   SUMMARY  OF  HEAT  TRANSFER  AND  FRICTIONAL  RESULTS: 
60x60  MESH  STAINLESS  STEEL  SCREEN  MATRIX 


HEAT  TRANSFER  RESTTT.TR. 


RUN 

1 
2 
3 
4 
5 
6 
7 
8 
9 
10 


REYNOLD1 
NUMBER 

431.57 

234.81 

159.26 

129.69 

106.78 

79.88 

50.69 

29.60 

15.43 

8.49 


NSTNPR 

.0267 

.0425 
.0521 
.0620 
.0659 
.0755 
.1040 
.1535 
.2584 
.4760 


CONDUCTION 
PARAMETER 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


MAXIMUM 
SLOPE 

.58333 

.68566 

.74639 

.80395 

.82926 

.87817 

1.05579 

1.22105 

1.57029 


N, 


tu 


0.0 


3.12 
4.96 
6.09 
7.25 
7.77 
8.83 
11.90 
17.94 
30.19 
2.11897   55.65 


ISOTHERMAL  FRICTION  RESULTS; 
RUN 


1 
2 
3 

.  4 
5 
6 
7 
8 
9 

10 


REYNOLD'S 
NUMBER 

437.93 

238.27 

161.60 

131.60 

108.36 

81.07 

51.44 

30.04 

15.66 
8.62 


FRICTION 
FACTOR 

.68548 

.82196 

.89263 

.97400 

1.06410 

1.23243 

1.57979 

2.31547 

4.47090 

9.50629 


TABLE  II   SUMMARY  OP  HEAT  TRANSFER  AND  FRICTIONAL  RESULTS: 
PLATE-FIN  REFERENCE  MATRIX 


HEAT  TRANSFER  RESULTS: 


RUN 

REYNOLD'S 
NUMBER ' 

*R 

NqmNPR2/5   CONDUCTION  MAXIMUM 
bi  r*           PARAMETER    SLOPE 

3         ^ 

Ntu 

'Ntu 

(A=o) 

1 

505.17 

.00277     .00128     .52583 

1.68- 

1.68 

2 

433.26 

.00429     o00150     .55712 

2.60 

2.60 

3 

387.19 

.00505     .00167     .58044 

3.06 

3.06 

4 

322.08 

.00634     .00201     .62383 

3.84 

■  3.84 

5 

248.92 

.00842     .00260     .69377 

5.10 

5.10 

6 

187.25 

.01162     .00346     .79434 

7.04 

7.04 

7 

125.74 

.01761     .00515     .94678 

10.67 

10.44 

8 

75.04 

.02889     .00862    1.16341 

17.51 

16.21 

9 

42.68 

.05699     .01516    1.43365 

34.54 

25.07 

10 

23.81 

.08403     .02717    1.43851 

ISOTHERMAL  FRICTION  RESULTS: 

RUN  REYNOLD'S  FRICTION 
#          NUMBER     FACTOR 
NR           f 

1  512.59     .02197 

2  439.62     .02485 

3  392.86     .02795 

4  326.79     .03304 

5  252.57     .04210 

6  189.99     .05433 

7  127.58     .08205 

8  76.14     .12974 

9  43.30     .23164 
10     24.17     .39637 

50.93 

25.22 

58 


TABLE  III 


SUMMARY  OF  HEAT  TRANSFER  AND  FRICTIONAL  RESULTS: 
10°  SKEW  STAINLESS  STEEL  MATRIX 


HEAT  TRANSFER  RESULTS: 


RUN 

REYNOLD'S 
NUMBER 

NR 

N  N  2/5 
WST1NPR 

3 

CONDUCTION 

PARAMETER 

X 

MAXIMUM 
SLOPE 

Ntu 

(A=0) 

1 

805.2? 

.00375  " 

.00081 

.46533 

1.32 

1.32 

2 

643.32 

.00440 

.00101 

. 51088 

1.55 

1.55 

3 

517.86 

.00695 

.00126 

.55129 

2.45 

2.45 

4 

402.58 

.01073 

.00162 

.62114 

3.78 

3.78 

5 

305.00 

.01488 

.00214 

.70167 

5.24 

5.24 

6 

200.72 

.02271 

.00324 

.83894 

8.00 

7.98 

7 

119.58 

.03642 

.00544 

1.02840 

12.83 

12.46 

8 

68.74 

.06376 

.00953 

1.27448 

22.46 

19.62 

9 

37.61 

.10788 

.01728 

1.44693 

38.00 

25.53 

ISOTHERMAL  FRICTION  FACTOR: 


RUN 

REYNOLD'S 

NUMBER 

FRICTION 

FACTOR 

f 

1 

817.04 

.02791 

2 

652.71 

.03405 

3 

525.43 

.04092 

4 

408.46 

.04951 

5 

309.46 

.06345 

6 

203.79 

.09412 

7 

121.41 

.15138 

8 

70.00 

.27472 

9 

38.16 

.47276 

TABLE  IV  SUMMARY  OF  HEAT  TRANSFER  AND  FRICTIONAL  RESULTS: 
20°  SKEW  STAINLESS  STEEL  MATRIX 


HEAT  TRANSFER  RESULTS; 


RUN 

REYNOLD'S 
NUMBER 

NR 

N  N  2fo 

1NST  PR 

CONDUCTION 
PARAMETER 

MAXIMUM 
SLOPE 

Ntu 

Ntu 
(A«0) 

1 

835.38 

.00394 

.00078 

.47359 

1.36 

1.36 

2 

728.65 

.00414 

.00090 

.48768 

1.43 

1.43 

3 

636.57 

.00475 

.00103 

.52394 

1.64 

1.64 

4 

535.10 

.00805 

.00122 

.56474 

2.78 

2.78 

5 

412.85 

.01182 

.00158 

.63657 

4.09 

4.09 

6 

312.06 

.01524 

.00209 

.70256 

5.26 

5.26 

7 

208.65 

.02167 

.00313 

.81581 

7.48 

7.48 

8 

126.47 

.03676 

.00517 

1.02339 

12.69 

12.33 

9 

71.21 

.06593 

.00917 

1.28940 

22.76 

20.09 

10 

39.79 

.09441 

.01646 

1.39790 

32.59 

23.79 

ISOTHERMAL  FRICTION  FACTOR: 


RUN 

i 

REYNOLD'S 
NUMBER 

NR 

FRICTION. 

FACTOR 

f 

l 

848.11 

.03197 

2 

739.39 

.03471 

3 

645.95 

.03908 

4 

542.97 

.04465 

5 

418,94 

.05431 

6 

316.65 

.06856 

7 

211.73 

.10255 

8 

128.33 

.16015 

9 

72.26 

.28784 

10 

40.40 

.  .48428 

TABLE  V  SUMMARY  OP  HEAT  TRANSFER  AND  FICTIONAL  RESULTS: 

200  SKEW  NICKEL  MATRIX  (PERFORATED  PLATES) 
HEAT  TRANSFER  RESTTT.TS, 


RUN 

i 

l 

2 

3 

4 

5 

6 

7 

8 

9 
10 
11 


REYNOLD' 
NUMBER 

966.07 
835.43 
738.42 
619.79 
484.37 
362.97 
239.79 
145.03 
82.26 
46.22 
23.50 


N  N  2/5 
1NST^PR 

d 

.00900 

.01146 

.01219 

.01418 

.01802 

.02497 

.03517 

.05199 

.07879 

.08147 

.09134 


CONDUCTION  MAXIMUM 
PARAMETER    SLOPE 
A 


.00337 

.00390 

.00441 

.00526 

.00674 

.00898 

.01359 

.02243 

.03951 

.07039 

.13843 


.56394 
.60157 
.61776 
.64762 
.71212 
.81813 
.95653 
,14669 
,39808 
42079 


2.72 
3.42 
3.70 
4.24 
5.41 
7.58 
11.29 
19.65 


1.50225 


tu 

(*=0) 

2.72 

3.46 

3.68 

4.28 

5.44 

7.54 

10.62 

15.70 

23.79 

24.60 

27.58 


ISOTHERMAL  FRICTTON  FACTOR: 


RUN 


1 
2 
3 
4 
5 
6 
1! 
8 
9 

10 

11 


REYNOLD'S 
NUMBER 

HR 
980.29 
847.72 
749.28 
628.91 
491.63 
368.31 
243.32 
147.16 
83.46 
46.89 
23.85 


FRICTION 

FACTOR 

f 

.03316 
.03709 
.04116 
.04686 
.05802 
.07208 
.10258 
.16737 
.29667 
.49889 
1.09090 


TABLE  VI   SUMMARY  OF  HEAT  TRANSFER  AND  FRICTIONAL  RESULTS: 
"PARALLEL  PLATE"  BRASS  MATRIX 


HEAT  TRANSFER  RESULTS: 


RUN 

REYNOLD'S 
NUMBER 

NR 

N  N  2/5 

WST  PR 

i 

CONDUCTION 
PARAMETER 

A 

MAXIMUM 
SLOPE 

Ntu 

a=o) 

1 

683.97 

.00376 

.00896 

.49410 

1.44 

1.46 

2 

551.23 

.00527 

.01112 

.54861 

2.27 

2.40 

3 

437.55 

.00748 

.01401 

.59598 

3.22 

3.36 

4 

340.15 

.01038 

.01802 

.66717 

4.47 

4.63 

5 

255.77 

.01335 

.02390 

.73472 

5.75 

5.63 

6 

168.67 

.02148 

.03634 

.87166 

9.25 

8.71 

7 

102.18 

.03483 

.05999 

.97562 

15.00 

11.15 

8 

57.06 

i  - 

.10761 

1.07854 

- 

- 

9 

33.37 

- 

.18369 

.98829 

- 

-  ' 

ISOTHERMAL  FRICTION  FACTOR: 


RUN 

REYNOLD'S 

NUMBER 

NR 

FRICTION 

FACTOR 

f 

1 

694.19 

.02846 

2 

559.46 

.03454 

3 

444.09 

.04292 

4 

345.24 

.05240 

5 

259.52 

.06907 

6 

171.19 

.10680 

7 

103.71 

.17239 

8 

57.90 

.32370 

9 

33.86 

.53085 

TABLE  VII   SUMMARY  OF  HEAT  TRANSFER  AND  FRICTIONAL  RESULTS: 

"PARALLEL  PLATE"  NICKEL  MATRIX  (PERFORATED  PLATES) 


HEAT  TRANSFER  RESULTS; 

NornN- 


tUN 

REYNOLD'S 

NUMBER 

NR 

1 

819.49 

2 

647.21 

3 

518.89 

4 

400.50 

5 

304.49 

6 

200.68 

7 

121.15 

8 

67.53 

9 

37.92 

2/3 


srPR 
3 

.00590 
.00810 
.00988 
.01458 
.01891 
.02659 
.04039 
.05502 
.04945 


CONDUCTION 
PARAMETER 

MAXIMUM 
SLOPE 

Ntu 

Q=o) 

.00393 

.54781 

2.32 

2.36 

.00497 

.58864 

3.18 

3.24 

•00620 

.62923 

3.89 

3.95 

.00804 

.73180 

5.77 

5.83 

.01058 

.81941 

7.61 

7.56 

.01602 

.95459 

11.31 

10.63 

.02658 

1.16142 

21.61 

16.15 

.04768 

1.34786 

- 

22.00 

.08478 

1.27951 

- 

19.70 

ISOTHERMAL  FRICTION  FACTOR: 
RUN  REYNOLD'S  FRICTION 


NUMBER 

NR 
831.56 
656.72 
526.51 
406.49 
309.04 
203.63 
122.96 
68.54 
38.48 


FACTOR 
f 

.02395 
.02955 
.03577 
.04439 
.05682 
.08725 
.14242 
.26776 
.45874 


8 
o 

.400 
445 
488 
529 
568 
603 
637 
697 
748 
791 
827 
857 
882 
903 
920 
934 
946 
956 
964 
978 
987 
992 
995 

000 
000 

000 

000 

000 
000 

229.= 

-           —                        —              —              J     J     J 

996 

6  96 

*r96 
296 
056 
9S6 
156 

SS6 

oe6 

1  16 

669 
999 
699 
0^9 
9^9 
961 

c 
c 

J 

o 

646 

677. 

703 

726 
745 
761 
775 
787 
798 
807 
815 
831 
843 
852 
859 
865 
870 
874 
877 
882 
886 
889 

893 
894 

895 

497 
522 

542 
560 
589 
615 
640 
662 
662 
6  99 

727 

749 
75£ 

767 
764 
600 
612 

6  3" 
844 

649 

956 

670 
873 
875 
67" 
879 

885 

667 
889 
891 

693 

6  94 

895 

.384 
.425 
.459 
.487 
•  51  1 
530 
544 
561 
571 
582 
592 
603 
615 
626 
637 
647 
658 
668 
678 
701 
722 
742 
760 
776 
792 
805 
817 
840 
860 
878 
893 
907 
921 
932 
942 
960 
975 
987 
997 
005 
013 
019 
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o 

o 
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.893 
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.973 
.996 
1.017 
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1.073 
1.105 
1.134 
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112-1 
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.922 
.959 
.998 
1.034 
1.068 

1.133 
1.193 
1.249 
1.302 
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1.0 
I.I 
1.2 
1.3 

1.4 
1.5 

1.6 
1.8 
2.0 
2.2 

2.4 
2.6 
2.8 
3.0 
3.2 

3.6 
3.8 
4.0 
4.5 
5-0 

6.5 
7.0 
7.5 
6.0 
9.0 
10.0 
11. 0 
12.0 
13.0 

14.0 
15.0 
16.0 
18.0 
20.0 
22.0 
24.0 
26.0 
28.0 
30.0 
32.0 
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APPENDIX  A. 
DATA  REDUCTION  RELATIONS 

The  following  summary  of  some  of  the  data  reduction  relations  is 
given  to  clarify  some  of  the  definitions  of  the  geometrical  parameters, 
to  show  how  they  influence  the  deduced  results  and  to  facilitate  an 
understand  as  to  the  exact  procedure  required  to  produce  from  raw  data 
a  point  on  the  j  vs  N  and  f  vs  N  curves. 

The  three  most  significant  geometric  parameters  that  afford  a  means 
of  comparison  between  matrices  are  so  defined  that  there  exists  a  unique 
relation  between  them.   Independent  measurement  of  any  two  establishes 
the  third  or  serves  as  a  check  on  its  independent  measurement.   The 
three  primary  surface  geometry  parameters  of  a  matrix  are: 

1)  HYDRAULIC  DIAMETER: 

4rH  =  4f  L  (A-l) 

2)  POROSITY:  Matrix  Flow 

Void  Volume      A  L       A 
p  =  Total  Volume   =   -c-   =  —  (A-2) 

fr       fr 

3)  AREA  COMPACTNESS:   Total  Heat  Transfer 

n  Surface  Area 

p  =    Total  Volume         =  A  A  (A-3) 

fr 

NOTE:   If  equation  (A-2)  is  divided  by  (A-3)  the  unique  relation  mentioned 
previously  is  formed. 

rH  =  p/0  (A-4) 

MAXIMUM  SLOPE: 

As  previously  discussed,  through  the  work  of  Locke,  Mondt,  Crestwick 
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.  rd,  it  has  been  established  that  the  maximum  slope  of  a  general- 
ized heating  or  cooling  curve  is  a  unique  function  of  N  and  /i  ,  the 
conduction  parameter, 


fcf2    "   fci 

d(t        -    t    } 
fl          i 

d(r/Ntu  )_ 

X  ) 


T  =  generalized  time  variable 


T  - 


t   =  initial  matrix  temp. 


N    =  0-±- 

tu    mc 


t   =  air  temp,  upstream 


air  temp,  downstream 


therefore : 


T/N_   =  mc  /W  c  Q 
tu      p   s  s 


and   d(  T  /N  )  =  mc  /  W  c   (d  0  ) 
tu      p    s  s 


.(A-5) 


also: 


and:  d(-~ ~ 

fl  "  Ci 

combining  (A-5)  and  (A- 6) 


d(  tf2  -  tfl) 


(A-6) 


d( 


") 


fl 
d(  T/N„  ) 


J 


d( 


•(A-7) 


The  matrix  is  heated  to  a  uniform  temperature,  t..   The  heat  source 
is  secured  and  the  difference  between  downstream  and  upstream  temperatures 
(t  ?  -  tf.)  is  recorded  vs.  time.   The  maximum  rate  of  change  of  this 
difference  with  time  may  be  determined  from  the  record. 

From  Figure  17: 


d(t 


12 


fl 


Chart  speed 


d© 

=  d&        sec 
y  =  d(t„  -  tpi) 


M 


Lf2 


Lfl 


fcfl' 


'l> 


inches 
inches 


Matrix  heat  Capacitance  _    TT     BTU/°F 
Fluid  heat  Capacitance       ~-   BTU/°F  sec 


Finally,  Equation  (A-7)  becomes: 


d(  T  /N,„) 


Chart  speed 


With  this  value  of  the  maximum  slope  and  a  calculated  value  of  A 
Howard's  curves  or  tables  may  be  entered  to  obtain  a  value  of  N 


MASS  FLOW  RATE: 

The  mass  flow  rate  is  calculated  as  outlined  in  reference  (7): 
m  =  359  K  Dq  FaY  |  A  PQ  t 
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where  K  =  flow  coefficient,  velocity  of  approach  factor 

Fa  =  thermal  expansion  factor 

Y  =  expansion  factor 

D  =  orifice  diameter,  inches 
o  * 

"^  =  specific  weight  of  fluid  flowing  at  the  orifice 

3 
temperature  lb/ ft  ,  corrected  to  existing 

atmospheric  pressure 

AP  =  Pressure  drop  across  theorifice  ,  inches  lUo 

REYNOLD'S  NUMBER: 

By  definition,  the  Reynolds  number  is 

4r  G 
NR=    ^—  (A-8) 

Where  G  is  the  mass  flow  velocity  based  on  the  flow  mean  cross- 
sectional  area,  A  . 


(A-9) 


c        P  fr 

Note,  if  (A-9)  is  substituted  into  (A-8) 


NR   = 

4m 
>"Afr 

but    from    (A-4) 

rH/n 

.  J- 

(  -7T  )  (A-10) 


-d) 


j3         "  R  juk 

Therefore,  for  a  given  range  of  mass  flows,  the  Reynolds  number  range 
for  a  given  matrix  is  proportional  to  the  inverse  of  its  compactness. 

nr  «:  Yfi 

n             A 
For  ease  of  calculation,  recall  from  (A-3)     H  -     — — 

fr 
substituting  into  (A-ll)  yields 
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(A-ll) 


4Lm 
A/4 


(A-12) 


THE  FANNING  FRICTION  FACTOR: 

The  following  sketch  describes  the  flow  system  to  be  considered. 


mmmznzzzzznnL 
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For  gas  flow  heat  exchanger  application,  the  pressure  changes 
from  section  1  to  a  and  from  b  to  2  are  very  small,  relative  to  the 
total  pressure,  and  therefore    "Vo,  ~  'Vj'     and  V^   ~  "l^. 

The  relation  for  the  flow  stream  pressure  drop  calculation  for 
most  heat  exchanger  cores  is  (8): 


P,    *jc  P, 


i*i-f)**fa-i)  +  f^-0-?-Vfi 


entrance       flow       core 
effect   acceleration   friction 


exit 
effect 


Solving  (A-13)  for  f  and  recalling  from  (A-l)   A/A  =  L/r  also 

I 
nT=  JQ         ,  we  find: 


■f- 


~^r~m-(fr*kh-Mcrfto+r. 


(A-13) 


(A-14) 


Now  the  perfect  gas  law  may  be  applied  since  testing  is  done  with  air  at 
moderate  temperatures  and  pressures.   Equation  (A-14)  becomes  for  the 
isothermal  case: 
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u 


where 


^Kt?)- ¥(^t)-(¥)(^o^) 


T  <A"15> 


J.  _  a        ,o   .  Pudfk 

(°     "       P  /     m  2 


and        T,    =  T„ 


The  first  term  is  by  far  the  largest  contributor  to  the  friction  factor 
for  small  pressure  differentials,  therefore  making  the  approximations 


ll   T   z2 


Equation(A-15)    becomes 


i-- 


'ifMf?)-(^vf  o+e1) 


(A-16) 


where : 


Kc  =  entrance  coefficient 
Ke  =  exit  coefficient 
whose  values  are  determined  from  appropriate  curves  in  Chapter  IV,  refer- 
ence (8). 

This  is  the  equation  used  to  determine  the  friction  factor  as 
presented  in  this  thesis. 

If  a  further  assumption  is  made  such  that 

{  *  Zy.  p.  ££  £  (a-17) 


and  equations  (A-4)  and  (A-9)  are  substituted,  Equation  (A-17)  becomes: 

1 
-JL  (A-18) 

(3 
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Therefore,  for  a  given  matrix  the  friction  factor  is  directly  pro- 
portional to  the  cube  of  the  porosity  and  indirectly  proportional  to  the 
compactness. 

f    *       P3//3  (A- 19) 

COLBURN   j    FACTOR: 

We  have  A    and  maximum  slope  with  which  we  get  N   from  Table  VIII 
or  Figure  18  taken  from  reference  (5). 

The  j  factor  is  defined  as: 

j  =  Norp  N   2/3  =  — N  n  2/3  (A-20) 

J    ST  PR         G  c     PR 
P 

substituting  (A-9)  for  G  and  multiplying  by   /A 

2/3  (A-21) 

'"WP 

by  definition: 

Therefore:  j  =   Ntu   A^   NPR 

Note,  from  (A-l)    Ac/A  =   rH/L,  substituting  f  /a  for 

r   from   (A-4)  gives: 

j  ■  Ntu  % 2/3  -r-    f'fi  (A-23) 

which  leads  to  the  relation  that  for  a  given  matrix,  j  is  directly  pro- 
portional to  the  porosity  and  inversely  proportional  to  the  compactness. 


hA                Ac 
mc                 A 
P 

"     NPR 

1        =       ?£- 

tu            mCp 

tu       A 

PR 

j  oc 


f'fi 


HEAT  TRANSFER  POWER  AND  FLOW  FRICTION  POWER: 

The  heat  transfer  power  per  unit  area  per  degree  temperature  differ- 
ence is  derived  as  follows: 
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J  ■  «„«.  2/3 


ST  PR 
hD, 


oM 


-L-       .      N    2/3 

NR 


Solving  for  h: 


>M  l 


2/3 
N„„  H 


ND  J 


Now  if  c  ,  ^/A     and  N   are  evaluated  at  standard  properties  which  by 
convention  are  taken  as  those  for  dry  air  at  500 °F  and  one  atmosphere; 

c      =  0.2477      BTU/lb    °F 
P 

yU  =  0.0678     lb/hr  ft 
/O    =  0.0413     lb/ft3 
N„„  =  0.671 


Thus,    the  heat   transfer  power  per  unit   area  per  degree   temperature 


BTU/hr   ft2°F 


The    flow   friction  power  per  unit   area    is   derived  as    follows: 
*1 


f    S     ilL^ft) 


1  uj     2>^ 
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E=AP(£)v(4>)--f 


v£>«\j 


where    6  =  p  V 


i(4f(f on«'] 


Now  ifyU  and   O  are  evaluated  at  standard  properties  (previously 
defined)  the  flow  friction  power  per  unit  area  evaluated  at  the  above 


ESTD  -   1.11  x  10-7  tjL)  3   (\o0>3   t   [hp/ft2J 

Figure  16  provides  a  comparison  of  the  20°  skew  and  the  "parallel 
plate"  matrices  from  a  heat  transfer  power  vs.  friction  power  point  of 
view  (h    vs.  EQTn).   For  comparison  purposes  the  surface  geometries 


STD 


STD 


«-3- 
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APPENDIX  B. 
DESCRIPTION  OF  EQUIPMENT 

To  obtain  heat  transfer  data,  a  transient  technique  was  employed. 
Briefly,  this  method  consists  of  heating  the  test  matrix,  to  a  uniform 
temperature  and  then  subjecting  the  matrix  to  a  step  change  in  fluid 
flow  temperature.   The  fluid  temperature  downstream  of  the  matrix  is 
recorded  versus  time. 

Friction  factor  was  calculated  from  pressure  drop  data  obtained 
from  static  pressure  taps  located  in  the  test  section  immediately  up- 
stream and  downstream  of  the  matrix. 

The  equipment  necessary  to  complete  the  experiment  falls  into  one  of 
the  following  systems. 

1)  Fluid  source 

2)  Flow  metering  system 

3)  Temperature  measuring  system 

4)  Pressure  measuring  system 

5)  Fluid  heater 

6)  Matrix  holder  &  test  section  casing 

FLUID  SOURCE: 

The  working  fluid  was,  of  course,  air.   Air  was  provided  to  the 
apparatus  by  placing  the  entire  rig  on  the  inlet  to  a  lOOcfm,  6  stage 
centrifugal  air  compressor.   The  compressor  requires  an  operating  voltage 
of  120/  at  90amp. 
FLOW  METERING  SYSTEM:  (See  Figure  1) 

The  flow  metering  system  consists  of  an  orifice  in  a  3.08"  inside 
diameter  aluminum  tube.   The  orifice  section  was  constructed  in  accord- 
ance with  the  ASME  specifications  of  reference  (7)  for  standard  D  and 
D/2  pressure  taps.   A  large  flow  range  was  obtained  by  means  of  7  orifice 
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plates  of  throat  diameter  2.310,  1.971,  1.232,  1.540,  .462  and  .308 
inches,  respectively.   Accurate  control  of  the  pressure  drop  across  the 
orifice  was  maintained  by  a  large  sliding  valve  downstream  of  the  ap- 
paratus on  the  inlet  to  the  compressor  and  another  sliding  valve  on  the 
discharge  side  of  the  compressor.   There  is  an  additional  valve  on  the 
inlet  side  of  the  compressor  that  is  fully  open  or  fully  closed.   It  is 
normally  opened  at  flow  rates  less  than  150  lb/hr  to  prevent  stalling  the 
compressor. 
TEMPERATURE  MEASURING  SYSTEM:  (See  Figure  1) 

Without  exception,  all  temperatures  were  measured  with  iron-constan- 
tan  thermocouples.   A  single  thermocouple  measured  the  orifice  tempera- 
ture (T  )  just  upstream  of  the  orifice  at  the  inlet  to  the  apparatus, 
o 

This  temperature  was  read  in  millivolts  on  a  RUBICON  CO.  portable  pre- 
cision potentiometer,  converted  to  temperature  and  recorded  for  each  run. 
Measuring  T   is  a  group  of  5  thermocouples  in  seriess  bound  together  and 
insulated  from  each  other  by  teflon  tape  and  inserted  in  a  1/8"  diameter 
aluminum  tube.   The  tube  was  placed  in  the  air  stream  with  a  small  window 
cut  in  it  facing  upstream.   At  this  point  ambient  air  temperature  was 
measured.   T«  and  T,  are  5  wire  thermocouple  grids  permanently  mounted  in 
the  test  section  cabinet.   Their  wires  are  connected  in  series  to  magnify 
the  thermocouple  output  voltage.   T_  is  a  movable  5  wire  thermocouple 
grid  located  in  the  matrix  holder  and  is  placed  adjacent  to  and  down- 
stream of  the  matrix,  it's  wires  are  also  in  series.   It's  position  is 
apparent  in  Figure  21.   The  thermocouples  arranged  in  this  manner  provided 
flexibility  and  complete  coverage  of  desired  points  of  interest.   If  the 
thermocouple  output  voltages  from  T~  and  T.  are  bucked  against  each  other, 
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the  temperature  measuring  system  response  curves  may  be  generated  on  a 
Leeds  &  Northrup  "Speedomax"  temperature  recorder.   T„  vs  T^  gives  the 
indication  of  temperature  uniformity  across  the  matrix  and  T.  vs  T 
records  the  rate  of  change  of  the  difference  between  upstream  and  down- 
stream temperatures  with  time.   It  is  this  record  that  provides  the  pri- 
mary data  for  the  evaluation  of  heat  transfer  data.   The  precise  relation 
of  this  record  to  the  Colburn  j  factor  is  explained  in  Appendix  A.   The 
record  is  made  on  a  Leeds  and  Northrup  "Speedomax"  Temperature  Recorder. 
The  recorder  has  variable  chart  speeds,  which  permits  maintenance  of  accept- 
able cooling  curves  with  regard  to  slope  determination.   The  Speedomax  also 
has  a  range  of  0-20  mv.   For  this  experiment,  the  recorder  is  precalibrated 
on  a  3  mv.  scale. 

There  is  one  addition  thermocouple  grid  that  is  not  used  in  the  experi- 
ment, but  was  used  in  checking  out  the  apparatus.   This  grid  would  take  the 
place  of  T  and  its  purpose  was  to  determine  the  temperature  distribution 
across  the   cross  section  in  both  the  vertical  and  horizontal  directions. 
It  was  placed  in  a  special  holder  that  also  held  a  velocity  probe  (See 
Figure  22)  and  was  lined  with  balsa  wood  to  provide  a  smooth  and  continu- 
ous internal  flow  passage. 
PRESSURE  MEASURING  SYSTEM:   (See  Figure  1) 

Pressure  taps  are  located  upstream  and  downstream  of  the  orifice  and 
upstream  and  downstream  of  the  matrix.   Each  tap  is  connected  by  Saran 
tubing  to  an  appropriate  manometer  or  draft  gage.   The  following  instru- 
ments were  used: 

1)  Type  "C"  micromanometer,  E.  Vernon  Hill  &  Co.  0-1.25" 

2)  Ellison  Differential  Direct  Draft  Gage,  0-6",  0-8" 

3)  Ellison  Inclined  Draft  Gage,  0-6" 

4)  Meriam  Instrument  Co.  0-60". 
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Any  one  or  combination  of  these  instruments  can  be  used  to  measure 
the  differential  pressures  of  the  orifice  and  the  matrix.   Frequent  cross 
checks  were  made  to  ensure  reliable  operation. 
FLUID  HEATER  SYSTEM:   (See  Figure  1) 

The  air  heater  section  consists  of  28  nichrome  .0031"  diameter  wire 
heaters.   The  heater  system  is  designed  to  heat  1000  pounds  of  air  per 
hour,  20°F  above  ambient  temperature.   Nichrome  was  chosen  because  of 
its  high  resistivity,  low  thermal  conductivity  and  low  specific  heat. 

The  nichrome  wire  is  wound  on  a  bakelite  frame,  two  heaters  to  a 
frame.   The  wires  are  1/32"  apart  and  there  are  50  and  52  wires  respective- 
ly, to  each  heater  on  a  frame.  The  heaters  on  each  frame  are  connected  in 
parallel  and  then  connected  through  a  switch,  to  a  variable  voltage  bus. 
There  are  two  frames  to  a  plastic  holder,  as  shown  in  Figures  23  and  24. 
In  these  figures,  the  heater  frame  and  its  electrical  connections  can  be 
clearly  seen. 

A  schematic  wiring  diagram  for  one  frame  is  shown  below.   All  14 
frames  are  wired  in  this  manner,  putting  all  28  heaters  in  parallel  with 
the  capability  of  switching  on  or  off  two  at  a  time. 
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The  28  heaters  arranged  in  this  way,  provide  a  variety  of  flexi- 
bility.  The  number  of  heaters  in  use,  should  decrease  as  the  flow  rate 
decreases.   Voltage  variation  can  be  obtained  by  a  GENERAL  RADIO  COMPANY 
"VARIAC",  0-240V,  8a,  50-60  cycles.   It  is  recommended  that  a  voltage  of 
220   or  in  this  vicinity,  be  maintained,  at  all  flow  rates,  and  the 
number  of  heaters  used,  decreased  until  the  number  of  heaters  is  four. 
At  this  point,  the  voltage  should  be  decreased,  as  necessary,  to  obtain 
a  20°  temperature  rise,  the  number  of  heaters  remaining  constant.   A 
note  of  caution:   Do  not  energize  a  heater,  unless  the  air  is  flowing, 
otherwise,  the  wires  will  sag  due  to  thermal  expansion  and  short  circuit- 
ing will  occur. 
MATRIX  HOLDER  AND  TEST  SECTION: 

The  matrix  holder  and  test  section  casing  are  made  from  polyethylene 
plastic.   This  material  was  chosen  because  of  its  low  heat  capacity. 

The  test  section  casing  holds  thermocouples  T«  and  T,  permanently 
mounted  upstream  and  downstream  respectively.   It  also  contains  the  up- 
stream and  downstream  static  pressure  taps  used  in  determining  the  fric- 
tion factor.   The  parts  making  up  the  casing  were  machined  to  close  toler- 
ances, to  ensure  a  tight  fitting  matrix  holder,  and  a  minimum  of  air  leaks. 

The  parts  comprising  the  matrix  holder,  as  implied  above,  were  machin- 
ed to  close  tolerances  to  ensure  a  tight  fit.   Great  care  was  taken  to 
ensure  smooth  alighment  of  the  flow  channel  through  the  heaters  and  the 
test  section.   The  size  of  the  matrix  holder  was  arbitrarily  chosen. 
The  flow  channel  dictated  the  width,  and  the  length  is  such  that  matrices 
of  various  flow  lengths,  up  to  3",  may  be  tested.   The  flow  channel  is 
3-3/16"  x  3-3/16". 

The  matrix  holder  has  one  special  provision  and  that  is  for  the  movable 
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thermocouple  T_.   A  small  hole  is  drilled  in  the  front,  and  the  wires 
are  led  through.   The  matrices  are  placed  in  the  holder  and  held  in 
place  by  bits  of  Styrofoam  plastic  cut  to  the  size  necessary  to  hold 
the  matrix  square  in  the  flow  channel. 

That  completes  a  description  of  the  various  systems  involved  in 
the  apparatus,  but  there  is  one  other  item  that  bears  mentioning.   Just 
upstream  of  the  heater  section,  is  a  section  provided  for  screens  to 
straighten  the  velocity  profile.   Provisions  were  made  for  1  to  4  screens. 
One  60  mesh  stainless  steel  screen  was  found  to  be  sufficient,  to  straight- 
en the  flow  at  the  flow  rates  used  (40-1000  lb/hr.) 
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APPENDIX  C 


CALCULATION  OF  GEOMETRIC  CONSTANTS 
FOR  PERFORATED  MATERIAL 

As  demonstrated  in  Appendix  A,  the  surface  geometry  has  a  strong 
influence  on  the  heat  transfer  and  friction  factor  curves.   Specifically, 
the  effects  of  p  and  ft   were  examined. 

In  using  perforated  materials,  the  presence  of  the  slots,  require 
modification  of  these  two  parameters.   It  is  the  intent  of  this  Appendix, 
to  demonstrate  the  calculation  procedure  utilized  to  determine  the  geo- 
metric parameters,  on  which  the  graphical  and  tabular  results  for  the 
matrices  constructed  from  perforated  material  are  based. 

The  material  used  in  the  perforated  plate  matrices,  was  a  pure  nickel, 
type  160/40TV,  slotted  opening  plate  material,  .0022  inches  thick.   This 
material  is  a  product  of  Perforated  Products,  Inc. 

The  basic  geometry  of  the  plate  is: 


/////Ik        /////A        /////A" 


[*  .036  cm  4*  .o24cmJ 


J  002} 


Using  the  average  between  the  top  and  bottom  of  the  tapered  slots, 
consider  an  average  slot  pattern  as: 

|  .023ovi  L .040  cm  _^> 

i 

2  cm 
1 
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CROSS  HATCHED  AREA  A 
c 

AREA  OF  SLOT  IN  A  .  = 
ch 

Now  define  y   =    Aslot 

ch 

V  _  'Q0QQ8 

0    ~      .001008 
Therefore,  plate  Porosity  (po)  =  (1 
Now  define 


(.063)  (.016)  =  .001008  cm 

2 

m 

=  AREA  REDUCTION  DUE  TO  SLOTS 


7.93% 


^      )        =    .9207 


r    _     CROSS    SECTION  AREA  OF   SLOTS 
"    "  CROSS    SECTION  AREA 


CONDUCTION  AREA 
REDUCTION  DUE  TO 
SLOTS 


r  _      (.040)(.0Q22)      _ 

°  (.063)(.0022)  J'D/° 

The  following  calculations  are  for  the  20°  skew  perforated  plate 


DIMENSIONS  OF  PLATE:  4.285"  x  1.98"  x  .0022" 

NUMBER  OF  PLATES-  142 

DIMENSIONS  OF  MATRIX:  3.255"  x  3.248"  x  1.98" 

WEIGHT  OF  MATRIX:  275.5735  grams 

MATERIAL  CONSTANTS;  A   =  555  lb/f  3 

V  38.7  BTU/hr.p  £t 

c   -    .1065   BTU/, 


'lb°F 

PLANE   SURFACE  AREA  =  (1«)  (4.  285)  (1.98)  (2)  =   ^33   ft2 

144 

TOTAL  HEAT  TRANSFER  AREA    (A)    =   PLANE    SURFACE  AREA   -   AREA  OF   SLOTS 
=    (1-  ^)      SURFACE  AREA 
=   (.9207)(16.733) 


A         =   15.41    ft 


FRONTAL  AREA    (A      )    =   3.255   x  3.248   =   10.572    in 


FREE  FLOW  AREA  =  FRONTAL  AREA  -  CROSS  SECTIONAL  AREA  CF  METAL 


=  .06416  ft 


CONDUCTION  AREA    (A   ) 


=    (142)(4.285)(.0022)(l-&) 

=    (142)(4.285)(.0022)(.375)    -    .50198   in2 


=      .003486   ft 


MAXTIX  VOLUME    (V   )    =  3.255"   x   3.248"   x   1.98"   =   20.933   in3 

m 

MATERIAL  VOLUME  (V  )  =  142x4. 285"xl. 98"x. 0022"  ^  2.650  in 

WEIGHT  OF  MATERIAL  (W  )  =  275.5735  gms./453.6  gm/lb  =  .6075  lb. 

MATERIAL  VOLUME  (\T)  =.  W  /  n   =  .6075  lb/555  lb/ft3  =  .001095  ft3 
S      s  A-*  s 

The  porosity  of  a  matrix  has  been  defined  as: 


MATRIX  FLOW  VOID  VOLUME 
MATRIX  VOLUME 


V   -  V   =  1  -  — 


p  =  1  - 


.01211 


.8733 


P  =  1 


.01211 


=  1  -  .09042  =  .90958 


The  compactness  of  a  matrix  has  been  defined  as 

O       TOTAL  HEAT  TRANSFER  AREA 
P  =  .       MATRiy  VOLUME 

Here,  we  will  redefine  A     to  include  the  effect  of  area  reduction: 


?- 


(PLANE  SURFACE  AREA) Q -  V  ) 
MATRIX  VOLUME 
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For   V  =  0 

P^^M^     -  "81.75 
For    )f  =  .0793 


5   _   (16.733)(.92Q7)  = 
P  "      .01211 


1272.18    ' ft 


The  hydraulic  diameter  is  conventionally  defined  as: 


.      FLOW  CROSS  SECTION  AREA 
H      ~    WETTED  PERIMETER 

Appendix  A  shows  that  hydraulic  diameter  can  also  be  expressed  as: 

"h  -  "  *'P 

To  be  consistent  with  the  conventional  definition,  for  the  case  of 
a  perforated  material,  A   is  used  vice  p       ,  since  the  slots  have  no 
effect  on  the  hydraulic  diameter. 

Finally  then, 

.8733 


D„  =  4 


.002528   ft 
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APPENDIX  D 
SAMPLE  CALCULATIONS 
lhe  following  sample  calculations  are  based  on  recorded  data  from 
Run  -v-3  for  20°  skew  matrix  constructed  from  perforated  nickel  sheets. 
This  is  chosen  primarily  becaust  the  basic  geometric  parameters  are  com- 
puted in  Appendix  C. 
RECORDED  DATA: 

D  =  2.310"  P    =  29.963  in.  Hg 

o  atm  & 

AP  =  2.61  in.  Ho0  CHART  SPEED  =  4.0  in/sec 

o  I 

W  s.6075  lb 

c  m 

H20 

T   =  68.0  °F 
o 

MATERIAL  CONSTANTS: 

c   =  .1065  BTU/lb  °F 
s 

k  =38.7  BTU/hr  ft°F 
s 

PRECALCULATED  GEOMETRIC   PARAMETERS:       (Re:      Appendix  C) 
A  =   15.41    ft2  p  =    .8733 


rJ1  "  ,o,o      ,o      ^2,^.3 


.0734   ft"  ft   =   1272.18    ft"/ftJ 


A     =   .06416   ft2  r     =   .000632    ft 

c  H 

A     =    .003486   ft2 
s 

L  =  1.98  inches  (  —  denotes  corrected  for  slot  effects) 
DETERMINATION  OF  MASS  FLOW:   as  outlined  in  Reference  (7) 

m  =  359  K  D    FY   \T~A~P   tf~" 
o   a     N     o 

where:    Y  Cz       1.0   Fig.  40b   Ref.  (7) 
F    «   1.0   Fig.  38   Ref.  (7) 


m    =      (359)(2.31)2  K    N  (2.6 


29.92 


=  849.40   K 
K  =  .749 


/S=      d/D 


2.31 
3.08 


.75 


m  =   (849.40)(.749)    =  636.21      lb/hr 
4m 
Rp      "  IT  VjU. 


Then 

Reynold's  number  for  the  pipe 


N„ 


(4)  (636.21)(12) 


=  7.35 
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"Rp  ~    (7T)  (3.08)  (.04410) 
From  Table  5,  Ref.  (7)  6=   0.75,   NR  =  7.35  x  104,  K  =  .749 


Therefore 


m  =  636.21  lb/hr 


MATRIX  REYNOLD'S  NUMBERS 

There  are  actually  two  Reynold's  numbers  to  be  calculated  for  the 
matrix.   One  is  for  isothermal  flow  and  is  associated  with  the  friction 
factor,  the  other  is  associated  with  heat  transfer.   They  differ  only  in 
the  temperature  at  which  the  absolute  viscosity  M.   is  evaluated.  yCA.  f 
(friction)  is  evaluated  at  ambient  temperature  (T  ).  JU.       (heat  transfer) 
is  evaluated  at  T  +  10°F  (an  estimate  of  the  average  bulk  fluid  tempera- 
ture) . 


N   )        =     — 
Vf  kfJi 


4Lm 

a>"h 


-    tO 


(4)(ffi   ft(636.21lbhr  _  61858 

(15. 41 ftz) (.04405  lb/ ft  hr) 

R'f       /*H  v ' ""*'    v    .04470' 


CONDUCTION  PARAMETER 


(38.7)(.003486)(12) 
(1.98)(636.21)(.24) 


.00536 


MAXIMUM  SLOPE 

From  Figure  26  the  slope  of  the  cooling  curve  was  determined  as 

described  in  Appendix  A.   The  slopes  of  three  curves  were  determined  and 

the  average  of  the  three  was  used  in  the  calculations. 


™™  -  r-  (SU    -  r- 


SLOPE  xCHART  SPEED 


C  =  m  c   =  (636.2llb/hr)(.24BTU/lb°F) 

=  .04241  BTI7sec  °f 
SLOPE  =  .16014  in"   (See  Figure  26) 
Finally: 


3600  sec 


MAXIMUM  SLOPE  = 


(.06470)   (.10775)  4) 


.65745 


(.04241) 
COLBURN  j  FACTOR 

With  A  and  MAXIMUM  SLOPE  enter  Table  VIII  and  find  N   =4.40. 

tu 

Linear  extrapolation  will  give  desired  accuracy. 
Now, 


J  =   N^„   A"  N 


2/3 


tu  A   PR 


from  Appendix  A,  Eq.  (A-22) 


15.41 

FANNING  FRICTION  FACTOR 

As  derived  in  Appendix  A  Equation  (A-16)  the  friction  factor  is: 


{- 


^§p)-(kc-O-£k(.  +  p0 


AP  =  3.96  in  Ho0  x  5.204  ?S *   - 
m  2  in  Ho0 


=  20.61  psf  =  .1431  psi 


psf 
in  H 

=  29.963  in  H  x  .4912  ?Sl„     =  14.7178  psi 


AF      =  4.69    in  Ho0  x  5.204      K „  n     =   24.41   psf   =    .1695   psi 
s  £  in   buu 


P  +P 

P     =    1,  2   =  14.4768  psi 
mean       2 


(P   „J  (144> 


^^eani_llZ:  =  (14.477)  (144)   =        lb   3 
fraean       RT  (53.3)(528)      -U^U/   ' rc 


.06416  ft2 
.000632 


H/L  "   1.98/12 


=  .00383 


K  =  -.738       K  =  .949   from  Chapter  IV  Reference  (8) 


O  APm  _   (64. 4)(. 07407)  (20.61)=   12.9579 

Tmean  2  7.587 


APm 


/i      2.  .  1431  psi  ,,    _,, . 

(1+P   >        "        14T4768Tsi      (1'761>    "    ' 


01742 


f  =  (12.9579  -  .01740  -  .211) (.00383)  =  (12. 7295) (.00383) 


.04859 


All  thermal  properties  (  U   ,   N  >k  ,  and  c  )  ,/ere  obtained  from 
References  (12)  and  (13).   The  temperature  at  which  the  properties  were 
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evaluated  was  the  ambient  temperature  with  the  exception  of  the  absolute 
viscosity  for  heat  transfer   (  /Liu    )• 
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APPENDIX  E 

A  FORTRAN  COMPUTER  PROGRAM 
FOR  DATA  REDUCTION 

A  short  computer  program  was  devised,  Figure  27,  so  that  from  basic 

experimental  data,  the  end  result  is: 

1)  Conduction  parameter  ( \) 

2)  Maximum  slope 

3)  Reynold's  number  for  heat  transfer 

4)  Reynold's  number  for  friction 

5)  Fanning  friction  factor 

It  was  anticipated  that  this  program  along  with  a  subroutine  would 
ultimately  produce  N   and  subsequently  the  Colburn  j  factor. 

Time  did  not  permit  the  adaption  of  Moreland's  work  (11)  to  a  sub- 
routine but  the  program  is  still  useful.   It's  value  when  data  for  a 
single  matrix  is  to  be  reduced  is  nebulous.   However,  the  reduction  of 
data  for  multiple  matrices  can  be  accomplished  with  a  considerable  savings 
of  time. 

The  following  is  a  list  of  inputs  necessary  in  Fortran  and  English 
language  along  with  the  required  dimensions: 
CARD  NO. 

30    XML      Matrix  length  ft 

2 
AS       Solid  conduction  area  ft 

RH       Hydraulic  radius  ft 

POR      Porosity 

2 
AFL      Flow  cross  section  area  ft 

2 
AHT      Matrix  heat  transfer  area  ft 

WM      Matrix  weight  lbs 
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CARD  NO. 


40    EK  Exit  flow  coefficient 

CK  Entrance  flow  coefficient 

SK  Solid  thermal  conductivity 

CM  Matrix  material  specific  heat 

N  Number  of  runs 

10    DO  Diameter  of  orifice  plate  opening 

DELPO  Pressure  drop  across  orifice 

GAM  Specific  weight  of  fluid 

ATMP  Atmospheric  pressure 

COR  Orifice  correction  K 

UH  Fluid  viscosity  for  heat  transfer 

UFR  Fluid  viscosity  for  friction 

10+2  CS  Chart  speed 

DELH  Pressure  drop  across  matrix 

HS  Static  pressure  matrix  inlet 

SLO  Slope  of  generalized  cooling  curve 

TEMPO  Orifice  temperature 
The  output  of  this  program  in  Fortran  and  English 

their  dimensions  is: 


BTU/hr°F   ft 
BTU/lb°F 


inches  H20 
lbs/ft3 
inches  Hg 

lb/hr  ft 

lb/hr  ft 

sec/inch 

inches  H_0 

inches  H-0 

(inches) 

°F 

language  along  with 


EMDOT  Mass  flow  rate  (m) 

RNUMP  Reynold's  number  for  the  pipe 

RNUMH  Reynold's  number  for  heat  transfer 

CAPFL  Thermal  capacitance  of  the  fluid  (C) 

CAP;  Thermal  capacitance  of  the  solid  (C  ) 

CONDPAR  Conduction  parameter  (  A  ) 


lb/hr 


BTU/sec  °F 
BTU/  °F 
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SLOMAX         Maximum  slope  of  generalized  cooling  curve 

2 
G  Mass  flow  velocity  lb/hr  ft 

RNUMFR         Reynold's  number  for  friction 

FFR  Fanning  Friction  factor 

Reduction  of  data  is  not  complete  since  j  remains  to  be  found. 

With  SLOMAX  and  CONDPAR,  enter  Table  VIII,  and  find  N   .   The  Colburn  i 

tu 

factor  follows  from  equation  (A-22),  Appendix  A. 

The  program  is  flexible  in  that  an  unlimited  number  of  matrices 
may  be  reduced  by  changing  the  first  card  of  the  program  appropriately. 
(L  =  the  number  of  matrices  to  be  reduced.) 
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